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x  SUMMARY 

■  ,•*. 

Extensive  measurements  have  been  made  on  a  plane  normal  to  the  wake  of  a  swept-wing/ 
fuselage  model  of  airbus  type,  at  =  0.49  in  incompressible  flow,  using  a  computer- 
controlled  wake-traverse  system  incorporating  a  null-reading  five-hole  yawmeter  probe. 

The  results  define  the  detailed  distributions  of  total  pressure  defect,  flow  velocity, 
flow  angles  and  streamwise  vorticity  in  the  viscous  wake  at  that  plane. 

The  results  have  been  analysed  within  the  basic  theoretical  framework  set  out  by 

Maskell  (1973),  allowing  the  calculation  of  lift  and  drag  from  traverse  data,  and  the 

resolution  of  the  drag  into  provisionally-defined  components  and  C'~>~  ,  respectively 

1  >  UII 

relating  to  the  profile  drag  and  vortex  drag. 
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I  INTRODUCTION 

The  work  described  in  this  Report  forms  part  of  a  long-term  study  of  the  wakes  of 
lifting  wings  and  wing-body  combinations  in  incompressible  flow.  It  is  based  on  measure¬ 
ments  made  with  a  computer-controlled  five-axis  wake-traverse  system,  employing  a  five- 

1  3-6 

hole  yawmeter  probe  in  a  self-aligning  or  null-reading  mode  ’  .  The  measuring  accuracy 

of  such  a  probe  used  in  the  conventional,  fixed  orientation  mode^  was  considered  insuf¬ 
ficient  in  the  present  context. 

The  primary  objective  of  the  current  work  is  the  trial  application  of  the  Maskell 
drag  analysis  scheme  to  a  comprehensive  set  of  total  pressure  defect,  flow  angle  and  flow 
velocity  measurements  on  a  plane  cutting  the  wake  of  a  model  representing  a  swept-wing 
airbus-type  transport.  Maskell*  reports  encouraging  results  from  preliminary  measurements 
in  the  wake  of  a  constant  chord,  uncambered  wing  (A  =  6)  on  a  simple  wire  support  rig. 

The  additional  complexities  of  camber,  sweep  and  taper,  fuselage  representation,  and  the 
use  of  a  three-strut  model  support  rig  are  now  introduced. 

I  .  2 

Maskell's  analysis  is  a  development  of  that  due  to  Betz  ,  and  allows  the  calcula¬ 
tion  of  the  lift  and  drag  of  a  wind-tunnel  model  from  measurements  confined  to  the 
viscous  wake  portion  of  a  wake-traverse  plane,  normal  to  the  tunnel  axis,  downstream  of 
the  model.  The  total  traverse-derived  drag  coefficient  (Cj^)  may,  subject  to  the  validity 
of  certain  assumptions  detailed  in  Ref  1,  be  resolved  into  two  components.  The  first 
(cDj)  is  broadly  consistent  with  the  conventional  concept  of  'profile  drag',  the  second 
(CDlI)  with  that  of  'vortex  drag'.  Maskell*  stresses  that  the  current  definitions  of 
Cpj  and  Cpjj  are  provisional  ones  only.  Their  physical  validity  and  their  uniqueness  - 
ie  their  invariance  with  the  chosen  location  of  the  traverse  plane  -  remain  to  be  fully 
substantiated  by  experiments  on  a  wide  range  of  model  and  support  rig  configurations. 

Ref  1  argues  the  potential  utility  of  this  drag  analysis  procedure,  once  its  sound¬ 
ness  and  its  generality  of  application  have  been  experimentally  established.  Once  one 
has  a  means  of  accurately  separating  the  profile  drag  and  vortex  drag  constituents  of 
lift-dependent  drag,  one  can  then  more  confidently  extrapolate  model-scale  test  results 
to  flight-scale,  particularly  for  high-lift  configurations.  It  will  then  be  possible  to 
re-examine  the  validity  of  the  planar  wake  assumptions  implicit  in  the  classical  vortex 
drag  prediction  theories,  and  to  study  three-dimensional  profile  drag  effects.  Again 
high-lift  configurations  will  be  of  primary  interest. 

A  secondary  objective  of  the  present  work  is  the  assessment  of  the  usefulness  of 
wake-traverse  measurements  for  diagnostic  purposes;  the  assumption  being  that  measurements 
would  only  be  made  in  a  limited  portion  of  the  viscous  wake,  and  made  in  less  detail  than 
that  required  for  drag  analysis.  It  is  reasonable  to  expect  that  such  measurements  could 
assist  in  the  identification  of  localised  sources  of  drag  that  might  be  considered  exces¬ 
sive,  such  as,  for  example:  wing-fuselage  junctions;  nacelles  and  external  stores;  flap 
and  slat  deployment  mechanisms;  wing  fences,  etc.  A  major  interest  in  the  present  work 
is  the  unfilleted  wing-fuselage  junction  of  the  model,  and  to  what  extent  the  junction 
vortices  (if  present)  contribute  to  the  profile  drag  and  vortex  drag  terms. 


Another  secondary  objective  of  the  present  work  is  the  assessment  of  the  extent  of 
mutual  interference  between  the  model  and  its  three-strut  support  rig.  The  wake-traverse 
procedure  clearly  cannot  yield  such  precise  quantitative  data  as  the  conventional  mechan¬ 
ical  balance  tares  and  interference  procedures.  It  can,  however,  be  expected  to  produce 
a  great  deal  of  useful  qualitative  information  as  to  the  mechanisms  determining  the 
viscous  component  of  the  mutual  interaction  between  the  model  and  its  support  rig. 

The  lift  condition  (Cj^  =  0.49)  in  the  present  investigation  is  comparable  to  that 
considered  in  Maskell's  pilot  study  of  the  wake  at  several  stations  downstream  of  a 

pitchwing*  ’*  though  the  current  Reynolds  number  is  some  4.5  times  higher  (Re-  =  950000) 

6  ^ 

The  current  wake-traverse  system  is  in  some  respects  more  advanced  than  Maskell's, 

particularly  with  regard  to  its  extensive  computerised  experiment-control  and  data- 

handling  facilities.  Its  greater  versatility  of  application  has  been  obtained  at  the 

cost  of  potential  increases  in  measurement  rate,  though  future  software  development  may 

recover  the  situation  somewhat. 

This  Report  also  discusses  some  aspects  of  the  results  of  some  preliminary  studies 
with  the  present  wake-traverse  system,  made  by  way  of  post-commissioning  trials.  Measure 
ments  were  made  at  three  lift  conditions  (one  corresponding  to  that  currently  considered) 
at  fixed  distance  aft  of  the  swept-wing  currently  considered.  The  only  configuration 
difference  was  the  absence  of  a  fuselage,  the  flaps  and  slats  being  retracted  in  all  four 
cases. 

Table  1  lists  the  configurations  studied  subsequent  to  Maskell's  small-scale  pilot 
study;  all  but  one  involved  cambered  swept-wings.  Tables  2  and  3  give  the  key  lift  and 
drag  results  for  the  three  wing-body  cases  analysed  in  detail,  while  Tables  4  and  5  give 
a  detailed  drag  breakdown  for  the  present  wing-body  case. 

Figs  1,  2  and  6a  show  the  present  experimental  configuration,  while  Figs  3  to  5 
show  the  basic  corrected  balance-derived  force  and  moment  characteristics  of  the  model. 

It  is  appropriate  to  note  that,  since  the  wake-traverse  drag  analysis  technique 
is  still  very  much  in  the  development  stage  and  low-speed  tunnel  manpower  very  limited, 
the  resources  have  not  been  available  to  effect  a  number  of  fairly  straightforward  com¬ 
puter  software  and  computerised  data  analysis  developments.  The  current  working  processe 
are  therefore  relatively  inefficient,  and  are  certainly  unrepresentative  of  a  fully- 
developed  wake-traverse  system  being  applied  to  a  long-term  drag  analysis  program.  At 
this  early  stage  it  has  been  felt  necessary  to  analyse  the  flowfield  data  in  very  con¬ 
siderable  detail,  rather  than  merely  submit  it  to  a  series  of  essentially  unmonitored 
computerised  integration  processes  with  the  sole  objective  of  deriving  the  values  of  the 
Cpj  and  CDjj  contributions  to  Cj3t  .  This  too  has  extended  the  timescale  of  each 
wake- traverse,  beyond  that  representative  of  routine  application. 

2  DETAILS  OF  THE  MODEL  AND  ITS  SUPPORT  RIG 

The  model  comprises  a  straight-tapered  28°  (quarter  chord)  swept-wing  of  aspect 
ratio  8.35.  The  wing  has  no  dihedral  or  twist  and  is  mounted  at  mid/low  position  on  a 
fuselage,  as  shown  in  Fig  1.  The  maximum  wing  thickness  is  10.7%  (at  37.5%  chord),  and 
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the  maximum  camber  is  1.155  (at  75Z  chord).  The  aerofoil  section  is  a  development  of  an 
NPL  type,  having  considerable  rear  loading.  Table  6  gives  the  basic  planform  dimensions 
and  Ref  8  gives  the  derailed  section  profile  co-ordinates. 

The  fuselage  is  representative  of  a  wide-bodied  transport  type,  and  the  wing  refer¬ 
ence  plane  is  at  1.1°  incidence  relative  to  the  fuselage  axis.  The  wing-fuselage  junc¬ 
tions  are  unlilleted,  and  the  diameter  of  the  circular  section  fuselage  remains  constant 
through  the  wing  root  region. 

The  model  has  been  tested  in  cruise  configuration,  with  the  slats  and  flaps  retrac¬ 
ted,  and  with  boundary  layer  transition  fixed  on  the  forward  fuselage  and  on  the  upper 
and  lower  wing  surfaces.  No  flap  hinge  or  flap  jack  fairings  have  been  represented  on 
the  wing,  and  all  joints  and  excrescences  were  carefully  smoothed  flush  with  the  surfaces. 

The  model  was  mounted  on  a  conventional  three-strut  support  rig  (Fig  2).  The  strut 
cleat  assemblies  on  the  undersurfaces  of  the  wings  and  of  the  rear  fuselage  were  not 
recessed  into  these  surface*. 

3  GENERAL  BACKGROUND  TO  THE  WAKE- TRAVERSE  ANALYSIS 
3. I  Theoretical  principles 

The  wake-traverse  drag  analysis  scheme  is  based  on  the  equation  for  the  conservation 
of  streamwise  (x)  momentum.  This  is  applied  to  a  control  volume  enclosing  the  model,  and 

otherwise  bounded  by  the  wake-traverse  plane,  the  displacement  surfaces  of  the  tunnel  wall 

boundary  layers,  and  a  plane  upstream  of  the  model  where  undisturbed,  basically  one¬ 
dimensional  flow  conditions  (HQ,  p^)  are  assumed  to  exist.  A  further  assumption  is  that 
no  streamwise  pressure  gradient  exists  in  the  empty  tunnel  working  section.  The  deriva¬ 
tion  of  the  model  lift,  using  the  same  control  volume  and  the  equation  for  the  conserva¬ 
tion  of  vertical  (z)  momentum  follows  essentially  similar  lines.  The  drag  and  lift 

equations  thus  derived  are  set  out  in  Fig  6b&c  respectively. 

The  lift  and  drag  forces  deduced  from  the  wake-traverse  measurements  relate  to  the 
constrained  flow  conditions  within  the  tunnel  working  section.  That  is,  these  forces  are 
in  principle  identical  to  those  sensed  by  the  tunnel  mechanical  balance  (where  one  is 
considering  the  model  and  the  wetted  surfaces  of  the  'live'  support  rig  as  an  entity); 
they  would  therefore  be  subjected  to  the  same  tunnel  corrections. 

The  basic  concept  of  a  scale-dependent  viscous  shear  layer  effect  ('profile  drag') 
and  an  independent  inviscid  secondary  flow  effect  ('vortex  drag')  is  a  useful  theoretical 
tool,  but  difficulties  arise  in  defining  these  effects  in  such  a  way  that  they  can  be 
isolated  in  terms  of  experimentally  measurable  quantities.  The  main  difficulty  is  that 
one  cannot  discriminate  between  that  portion  of  the  total  head  decrement  measured  on  the 
traverse  plane  that  is  attributable  to  'profile  drag'  (Cj) ^ )  effects,  and  that  portion 
attributable  to  the  viscous  dissipation  of  strearwise  shed  vorticity  between  the  wing 
and  the  traverse  plane,  ie  'vortex  drag'  (Cp.^)  effects.  For  example,  streamwise  vor¬ 
ticity  is  shed  within  the  wake  flow  inboard  of  the  tip  vortex.  In  this  region  it  is 
reasonable  to  assume  that  the  viscous  dissipation  of  streamwise  vorticity  will  be  neglig¬ 
ible,  the  velocity  gradients  typically  being  relatively  low.  This  is  not  the  case  for 
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the  tip  vortices,  however,  where  extremely  high  radial  velocity  gradients  exist  and  sur¬ 
face  shear  layer  fluid  from  the  wingtip  regions  is  entrained  into  the  vortices.  For  the 
purposes  of  the  analysis  o'  the  results  of  the  present  investigation  it  has  been  assumed 
that  the  total  head  defects  in  the  vortex  cores  (measured  as  a  'profile  drag'  contribution 
ACq.  )  are  dominantly  produced  by  the  dissipation  ot  (Ip  related  rotational  energy 

generated  in  an  idealised  inviscid  secondary  flow.  Thus  the  increment  £CD  has  been 

1  v  c 

accounted  to  the  Cp^  term. 

If  the  rate  of  viscous  dissipation  in  the  vortex  core  is  significant,  and  the  above 
4CpT  transfer  is  not  made,  then  one  would  expecL  in  interchange  between  the  measured 
Cp^  and  components  of  a  constant  Cp^  -  as  ..ie  traverse  plane  is  moved  progres¬ 

sively  further  aft  of  a  wing  at  given  lift.  The  Cp  term  would  incr  ase  with  progres- 

1 5  ‘ 

sive  dissipation  of  vorticity.  Maskell  found  just  such  an  effect  in  thL  analysis  of 
the  first  two  of  the  three  traverses  behind  a  pitchwing  (Ref  1  reports  the  first  traverse 
and  the  basic  theory),  but  was  unable  to  account  for  it  fully  with  the  above  ACp 

I  vc 

adjustment. 

Throughout  the  remainder  of  this  Report  the  terms  'profile  drag'  and  'vortex  drag' 
are  frequently  used  with  reference  to  the  measured  Cp^  and  cDjj  components.  The 
single  wake-traverse  dataset  considered  cannct  be  taken  either  to  confirm  or  to  disprove 
the  reservations  expressed  in  section  I,  as  to  the  validity  of  the  Cp^  and  ’-Du 
definitions. 


3 . 2  Practical  application  of  the  drag  analysis  procedure 


Fig  6a  shows  the  regions  of  the  wake-traverse  plane  within  which  measurements  have 
been  made.  The  arbitrarily  defined  analysis  regions  1-10  indicate  the  limits  of  measure¬ 
ments  on  rectangular  meshes  of  points.  Each  analysis  region  extends  a  little  way  beyond 
the  extremities  of  the  viscous  wake,  and  has  only  a  limited  variation  of  mesh  size  within 
it. 


Fig  6b  summarises  the  key  expressions  occurring  in  the  drag  analysis  and  shows  the 
quantities  (subscript  '2')  that  must  be  measured  at  each  mesh  point  on  the  traverse  plane. 
The  tunnel  reference  parameters  Up,  Hp,  qp  are  defined  at  the  centre  of  the  upstream 
reference  plane,  and  the  theory  assumes  these  to  represent  the  whole  upstream  plane.  This 
assumption  is  not  entirely  valid  (Fig  14)  but  appropriate  corrections  can  readily  be  made 
to  the  Cqj  or  'profile  drag'  term  (section  4.6). 

The  auxiliary  terms  ACq^  and  ACq^  (Fig  6b;  equations  (6),  (7))  represent 

induced  upwash  contributions.  One  can  only  measure  the  quantity  /T'^ds  on  the  traverse 

1  W 

plane,  not  the  quantity  /¥£ ds  employed  in  the  basic  flow  model.  T  is  the  stream 
function  of  the  transverse  flow  that  would  exist  far  downstream  in  wind-tunnel  if  the  dis¬ 
tribution  of  streamwise  vorticity  (£(y,z))  measured  in  the  wake-traverse  plane  were  to 
persist  unchanged  to  infinity  downstream.  ¥’  is  the  actual  stream  function  measured  in 
the  traverse  plane  for  the  approximately  two-dimensional  flow  in  the  near  field  of  the 
wing;  and  ACq^^  relate  T'  to  f  defines  a  direct  effect  dependent 

on  the  spanwise  loading  on  the  wing,  and  ACq  defines  a  wall  constraint  or  image 

1  lb 

effect.  In  evaluating  these  auxiliary  terms  x  has  been  measured  relative  to  the  mean 
quarter  chord  station  on  the  wing. 
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The  first  auxiliary  term  requires  a  circulation  distribution  T(y)  in  its  evalu¬ 
ation;  this  is  equivalent  to  the  surface  integral  of  streamwise  vorticity  in  the  viscous 
wake  on  the  traverse  plane,  outboard  of  each  of  a  series  of  stations  at  various  y 
values.  T ( y )  only  emerges  from  the  lift  calculations  discussed  below,  and  so  the  evalu¬ 
ation  of  the  auxiliary  term  cannot  be  incorporated  readily  into  the  computer  data  reduc¬ 
tion  scheme.  Both  terms  have  been  calculated  manually  and  appear  in  Table  4. 

A  correction  term  is  applied  to  the  term  (Fig  6b,  equation  (8)),  correc¬ 

ting  the  C[)j  overstatement  otherwise  produced  by  the  blockage  velocity  effect 

(Al'  =  U  -  U.  =  u,_)  of  the  presence  of  the  model. 

u  e  0  0 

3 . 3  Practical  application  of  the  lift  analysis  procedure 

Fig  6c  shows  the  key  expressions  involved  in  the  lift  analysis,  together  with  a 
schematic  representation  of  the  starboard  half  of  the  symmetrical  'liftloop'  employed  in 
tb°  calculations.  This  contour  lies  entirely  in  the  free-stream,  other  than  for  unavoid¬ 
able  crossings  of  support  rig  wakes  (segment  'L-K*  for  example);  for  convenience  it 
largely  follows  the  edges  of  the  arbitrarily  chosen  analysis  regions. 


The  dominant  contribution  to  the  local  lift  coefficient  (Fig  7)  at  a  spanwise 

station  'y'  is  pl^F  /(q^bc)  (Fig  6c,  equation  (2)),  where  T^  is  the  circulation  round 
the  contour  formed  by  the  line  cutting  the  wake  at  'y'  (segment  'A-P'  for  example)  and 
that  portion  of  the  liftloop  lying  to  its  right-hand-side.  T  is  calculated  thus, 
progressing  from  the  starboard  tip  to  the  port,  sequentially  treating  each  'y'  station 
for  which  measured  data  are  available.  The  liftloop  or  contour  integration  approach  is, 
for  the  present  purposes,  more  convenient  than  evaluating  the  surface  integrals  (for  the 
streamwise  vorticity)  for  the  areas  enclosed  within  the  set  of  contours  corresponding  to 
the  various  ' y'  values. 


It  has  been  found  that  the  use  of  rather  than  U  makes  a  negligible  differ- 

0  e 

ence  to  the  values  derived  from  this  procedure,  so  has  been  used  with  a  con¬ 

siderable  resulting  simplification  of  the  data  reduction  scheme. 

* 

The  contribution  from  the  second  term,  the  integral  P/(U-  -  U_)W_dx  along  the 

W  Z  L  J. 

cut  through  the  wake,  has  consistently  been  found  to  be  negligibly  small. 

The  third  term  in  the  equation,  namely  pU^AT^/  (^q^c)  corrects  for  the  circula¬ 
tion  jump  encountered  in  crossing  the  main  strut  wakes  (segment  'L-K'  for  example),  when 
calculating  in  the  region  -y'  <y  <  +y '  between  the  struts.  These  struts  are  in 

fact  lifting  surfaces,  given  the  sidewash  field  generated  by  the  wing.  This  factor  did 
not  require  consideration  in  Maskell's  application  of  the  theory  to  a  wing  on  a  wire 
support  rig'.  The  circulation  is  evaluated  taking  the  anticlockwise  direction  as 

positive,  and  the  circulation  increments  produced  by  the  streamwise  vorticity  shed  from 
the  wing  on  the  segment  'PONMLKJH'  (Fig  6c)  of  the  contour  in  the  strut  region  are  in 
fact  positive,  taking  the  starboard  strut  as  example.  The  circulation  jump  produced  by 
the  strut  wake  will  be  negative  in  terms  of  the  above  convention,  so  the  correction 
AT.  requires  to  be  positive.  The  same  argument  holds  for  the  port  strut  wake,  all 
signs  being  reversed. 
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The  strut  circulation  correction  AT  (=  +y  ,  Fig  8b)  has  been  evaluated  on  a 

L.  K  Z 

rectangular  contour  close  to  the  starboard  main  strut  at  station  'A'  (z  =  -334  isd)  ,  about 
half  way  along  the  wetted  strut  length.  The  streamlines  crossing  the  rear  segment  of  the 
contour  at  'A'  reach  the  main  wake  traverse  plane  at  a  height  z  =“  341.5  mm  ,  at  which 
point  the  wake  crossing  segment  L-K  (Fig  6c)  of  the  liftloop  has  been  located,  after 
having  made  the  necessary  supplementary  measurements.  (The  streamlines  were  traced  aft 
of  'A'  to  the  traverse  plane,  by  searching  for  the  local  high  total  head  defect  region 
generated  by  an  excrescence  temporarily  attached  to  the  strut  trailing  edge  in  plane  'A'.) 

It  has  been  assumed  that  the  port  strut  wake  has  an  identical  y^  to  the  starboard 
wake,  at  the  same  z  ,  a  reasonable  assumption  given  the  approximate  constancy  of  the 
absolute  circulation  ^CLg- )  noted  on  the  starboard  strut  in  the  region  of  'A'.  Accord¬ 
ingly  the  lift  correction  increment  Cl<,c  =  pU^Ar^j,/ (q^bc)  has  been  added  to  all  C^, 
values  computed  for  -y'  <  y  <  +y'  ,  as  shown  in  Fig  7.  A  smooth  curve  has  been  faired 
through  the  strut  regions  on  the  C^y  plot,  and  the  overall  traverse-derived 
obtained  by  spanwise  integration. 

It  should  be  noted  that  C^y  represents  the  spanwise  distribution  of  the  defect  of 
vertical  momentum  component  in  the  traverse  plane.  As  a  consequence  of  the  sidewash 
effects  generated  by  the  wing  in  its  wake,  and  also  of  flow  convergence  over  the  rear 
fuselage,  C^y  cannot  be  regarded  as  representing  the  spanwise  lift  distribution  actually 
existing  on  the  wing  (even  though  it  has  been  taken  as  such  for  the  estimation  of  some 
small  cDjj-  auxiliary  terms).  These  crossflow  effects  will  be  particularly  strong  in 
the  centre  section,  so  the  CLy  diagram  (Fig  7)  gives  no  directly  useful  information  on 
the  wing-body  interaction  effect  on  spanwise  loading. 

3 . 4  Measurement  of  the  strut  circulation  effect 

Measurements  have  been  made  close  to  the  starboard  main  strut,  on  contours  enclosing 
the  strut  section,  for  each  of  a  number  of  constant  z  planes,  using  the  wake-traverse 
gear.  Although  it  was  not  possible  to  make  flow  measurements  immediately  upstream  of  the 
strut  nose,  it  was  found  that  one  could  make  extrapolations  with  sufficient  confidence  to 
close  the  circulation  loop  and  thus  compute  y^  .  Fig  8b  shows  the  variation  of  crossflow 
induced  circulation  along  the  strut  thus  derived.  When  y^  is  expressed  in  f°rrn 

(ie  nondimensionalised  with  respect  to  local  strut  chord  C(Z)  )  the  sectional  circulation 
is  constant  over  a  considerable  portion  of  the  exposed  strut  length.  The  rapid  drops  in 
sectional  y near  the  lower  and  upper  ends  of  the  strut  can  respectively  be  attributed 
to  strut  guard  leakage  plume  effects  and  to  model  wake  effects. 

In  addition  to  establishing  y^  for  the  wing  plus  body  case,  the  exercise  included 
the  measurement  of  y  at  station  'A'  for  three  model  lift  conditions  for  the  wing-alone 

1 3  f 

configuration  .  This  data  (Fig  8a)  shows  that  the  strut  circulation  y^  ,  one  of  the 
interference  effects  of  the  model  on  the  support  rig,  is  strongly  dependent  on  the  lift 
developed  by  the  model.  At  C^c  a  0.49  the  strut  circulation  for  the  station  'A'  is 
some  50%  greater  for  the  wing  and  fuselage  case  than  for  the  wing  alone  case,  a  disparity 
which  is  far  greater  than  proportionate  to  the  probable  changes  in  distribution  on 

the  wing,  in  proximity  to  the  strut,  produced  by  introduction  of  a  fuselage. 


EXPERIMENTAL  DETAILS 


4.1  General  comments 

The  tests  were  made  in  the  No.  2  I  1}  it  *  81  It  (3.5  m  *  2.6  m)  low— speed  tunnel  at 
Earnborough,  at  a  windspeed  of  60  m/s  corresponding  to  a  mean  chord  Reynolds  number 
Re_  =  0.95  *  10h  .  The  model  was  tested  without  a  tailplane,  fin  or  engine  nacelles. 

boundary  layer  transition  was  f !  xed  at  the  5%  chord  line  over  the  full  span  of  the 
upper  wing  surface,  similarly  at  the  15%  chordline  on  the  lower  surface,  and  at  the  50% 
maximum  fuselage  diameter  station  on  the  model  nose  -  using  wires  in  each  case.  Some 
preliminary  flow  visualisation  tests  were  made  -  employing  an  acenapthe 1 ene/pe t rol eum  - 
ether  surface  deposit  -  to  confirm  tin-  effectiveness  of  the  transition  fixing,  particu¬ 
larly  at  the  chosen  wake-traverse  incidence  ( ■«  =5.10).  It  was  considered  essential  to 

precLude  the  possibility  of  transition  'creep'  over  the  considerable  duration  of  traverse 
measurements . 

4 . 2  Experimental  programme 

This  required  some  125  tunnel  running  hours  over  a  one  month  period.  It  comprised 

(i)  an  initial  set  of  mechanical  balance  measurements,  to  establish  the  basic 
characteristics  of  the  model  and  the  effect  of  transition  fixing  (FiRs  3  to  5), 

(ii)  a  detailed  traverse  of  the  viscous  wake  in  a  plane  7.4  tip  chordlengths  aft 
of  the  tip  t rai 1 i ng-edges  (Figs  I  and  6a),  measurements  being  made  at  5100  points, 
and , 

(iii)  several  partial  traverses  close  to  the  fuselage  side,  seeking  evidence  of 
discrete  vortices  or  excessive  total  head  losses  -  generated  in  the  wing/ fusel  age 
junction  region. 

The  force  measurements  were  made  over  a  range  of  uncorrected  incidence 
-5°  <  n.j,  <  1 2°  ,  the  lower  limit  being  fixed  by  separation  on  the  wing  under-surface 
(causing  model  buffeting),  and  the  upper  limit  by  stalling.  The  wake  traverse  measure¬ 
ments  were  carried  out  at  a  mean  corrected  lift  coefficient  Cj^  =  0.49,  corresponding 
to  a  corrected  incidence  =  5.25°(a^,  =  5.1)  .  During  traversing  the  test  condition 
was  monitored  by  periodic  mechanical  balance  measurements.  No  consistent  trends  of 
change  in  lift  or  drag  were  discernible,  within  the  t0.8%  scatter  band  of  lift  data  or 
the  H . 2%  band  of  drag  data  (relative  to  mean  values).  Concurrent  checks  on  the  nominal 

■> 

dynamic  pressure  produced  values  lying  well  within  *0.2%  of  that  desired  (q(1  =  2190.7  N/m~). 

4 . 3  Reduction  of  forces  data 

The  balance  measurements  have  been  reduced  to  conventional  coefficient  form,  using 

q 

a  standard  computer  program  .  For  the  purposes  ot  Figs  3  to  5  these  coefficients  have 
been  corrected  -  for  stream  pitch,  wall  constraint,  s'  Teamwise  pressure  gradient,  solid- 
blockage  and  strut  deflection  (model  CG  offset)  effects.  No  corrections  have  been  applied 
in  respect  of  separated  flow  wake-blockage  effects  (assumed  negligible),  or  of  the  mutual 
interference  between  the  model  and  its  support  struts'®. 
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4 . 4  Make— traversing  -  experimental  technique 

The  measurements  were  made  using  a  five-tube  yawmeter  with  a  central  total  pressure 
tube.  It  was  of  conventional  design'*,  constructed  of  0. 5  mm  diameter  (0.27  mm  bore) 
hypodermic  tubing,  and  having  a  nose  apex  angle  of  60°.  It  was  employed  in  a  null-reading 
mode,  being  aligned  with  the  local  velocity  vector  before  taking  each  set  of  reading  - 
namely  five  pressures  [  H,  Pj  ,  p^ ,  ,  p^l  and  the  probe  orient?*  .n,  angles  [  6 ,  t ;  relative  to 

the  tunnel  centreline]  -  at  a  given  probe  tip  location  I  x,y,z;  relative  to  the  intersec¬ 
tion  of  the  line  joining  the  trailing-edge  tips  and  the  plane  of  model  symmetry  (y  =  0)]  . 
Probe  alignment  was  effected  by  rotation  of  the  probe  support  quadrant  about  a  roll  axis 
(-200°  <j>  < +200°)  parallel  to  the  tunneL  centreline,  and  by  circumferential  movements 

of  the  circular  arc  quadrant  in  a  range  -30°  <  0  <  +30°  (or  alternatively  0°  <  9  <  60°) 
where  9=0°  corresponds  to  the  aerodynamic  axis  of  the  probe  being  parallel  to  the 
tunnel  axis.  The  6,<(>  (quadrant  and  roll)  motions  have  a  virtual  centre  at  the  probe 
tip. 

Lovell^  gives  a  detailed  description  of  the  mechanical  operation  of  the  traverse 
gear  and  the  system  providing  probe  nulling;  he  also  gives  details  of  the  alternative 
manual,  computer-assisted  manual,  and  fully— computer-controlled  inodes  of  operation  of  the 
wake- traverse  system.  The  system  allows  flow  measurements  to  be  made  at  seauences  of 
points  to-ming  prescribed  lines  or  grids  of  points  in  any  chosen  portion  of  the  tunnel 
working  section,  entirely  under  computer  control  where  appropriate. 

The  probe  nulling  motions  of  the  present  wake-traverse  system  differ  from  those 
described  by  Maskell*;  the  latter  corresponded  to  direct  movements  in  the  pitch  and  yaw 
senses,  rather  than  in  quadrant  and  roll  senses.  The  current  system  inherently  has  a 
somewhat  slower  rate  of  probe  nulling  -  and  therefore  a  slower  rate  of  measurement  -  for 
a  given  standard  of  control  system.  It  does  however  have  the  advantage  of  reduced  aero¬ 
dynamic  interference  on  the  flow  under  study  -  due  to  the  elimination  of  significant 
probe  support  asymmetries. 

The  complexities  of  the  mechanism  and  the  control  system  required  to  effect  nulling 
of  the  probe  can  be  eliminated  if  it  is  used  in  fixed  orientation.  One  then  has  to 
establish  calibrations  for  the  flow  angles  between  the  probe  axis  and  the  local  velocity 
vector,  and  the  total  head  coefficient  and  static  pressure  coefficient,  as  functions  of 
the  five  measured  pressures;  Ref  7  describes  such  an  approach.  However,  such  calibration 
schemes  have  been  claimed**  to  give  typical  accuracies  of  dynamic  pressure  measurement 
of  only  3-5£,  for  flow  angles  less  than  25°.  Such  an  accuracy  would  be  unacceptably  low 
in  the  present  drag  analysis  context. 

Fig  I  shows  the  location  of  the  main  wake-traverse  plane  relative  to  the  model,  and 
also  those  of  the  supplementary  planes  'A'  and  'B'  relative  to  the  fuselage  sides.  Fig  6a 
shows  the  division  of  the  traverse  plane  into  rectangular  analysis  regions,  each  extending 
a  short  distance  beyond  the  wake  extremities.  Each  analysis  region  is  composed  of  regular 
rectangular  grids  of  measurement  points,  chosen  according  as  to  the  local  gradients  of 
flow  angle  and  total  head  encountered,  as  follows. 
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Region  No. 

Region 

Ay  (mm) 

Az  (mm) 

1  .  3 

i  vortex  cores 

4 

4 

‘  vortex  periphery 

8 

8 

{  inner  wing 

40 

8 

2 ,  A 

\  strut  top/raid  wing 

8 

8 

f outer  wing 

40 

8 

3 

centresec t ion 

20 

8 

6,7,8 

struts 

8 

40 

Experience  has  shown  that  it  is  often  important  to  avoid  changes  of  grid  pitch 
within  a  reference  area;  small  economies  in  tunnel  running  time  can  be  offset  by  the 
introduction  of  difficulties  in  interpreting  the  local  vorticity  data  at  the  intersections 
of  meshes  of  differing  pitch.  Constant  values  of  grid  Ay,  A?.  are  particularly  useful  in 
the  centresec tion  reference  area  -  where  the  wingroot  vortices  may  have  migrated  further 
towards  the  centreline  than  expected  -  and  also  in  the  tip  vortex  reference  area,  where 
one  wishes  to  examine  the  manner  in  which  the  viscous  wake  is  rolled  round  the  vortex 
core.  In  the  present  tests  the  relaxation  of  the  mesh  size  on  the  vortex  periphery  -  on 
ground  of  economy  of  measurement  points  in  the  light  of  the  relatively  lower  gradients 
of  H,9,$  observed  -  has  proved  to  be  a  mistake.  It  led  to  a  loss  of  detail  in  the 
measured  vorticity  distribution  between  the  vortex  core  and  the  viscous  wake  rolled 
round  it.  The  consequent  loss  in  the  accuracy  of  measurement  is  negligible,  but 

one  has  lost  some  useful  data  on  the  wake  roll-up  process. 

During  the  wake- t rave rse  test,  the  tunnel  mechanical  balance  readings  (0^,0^)  and 
the  tunnel  speed  were  routinely  checked,  the  model  incidence  (aT  =  5.1°)  being  held 
constant  throughout.  Additionally,  measurements  were  made  at  a  check  point  before  and 
after  each  measurement  run  -  to  highlight  any  changes  in  test  condition  or  probe  calibra¬ 
tion  characteristics,  or  system  failure  during  the  run.  The  check  point  disparities  fell 
entirely  in  the  third  category,  and  the  relevant  data  has  been  discarded  where  necessary. 

The  check  point  was  located  in  the  free-stream  outboard  of  the  starboard  tip  vortex 

(where  the  local  yaw  angle  was  small),  and  it  was  also  used  to  set  the  zero  for  quadrant 

angle  (0)  measurements  -  by  'balancing'  or  nulling  the  probe  with  the  local  velocity 

o  •  •  o 

vector  at  two  roll  angles  approximately  180  apart.  The  deviation  (A$  -  180  )  is  a 

D 

function  of  the  inevitable  small  probe  asymmetries  and  also  of  the  ljcal  quadrant  angle 
(9fi),  where 


Because  of  the  existence  of  the  roll  asymmetry  the  'unstable'  balance  point  4’;^) 

can  only  be  located  by  reversal  of  the  usual  convention  that  positive  roll  error  signals 

(Ap  )  require  negative  increments  of  roll  angle  $  to  null  them.  The  roll  asymmetry  is, 
K 

for  the  purposes  of  probe  calibration,  characterised  by  the  roll  tube  angular  position 
error 
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£ 

P 


(2) 


where  6^  is  the  measured  quadrant  angle  (noting  that  one  sets  a  zero  with 

9  =  6g  =  -0b,  c  )  and  A<£R  the  roll  angle  increment  between  the  normal  (or  'stable') 

B  S  US  ® 

and  the  'unstable'  balance  points.  Following  the  recommendation  of  Ref  6,  e  has  in 

P 

fact  been  computed  as  the  mean  of  values  derived  for  a  number  of  6  values  correspond- 

D 

ing  to  a  range  of  points  in  the  flowfield. 

The  roll  angle  zero  is  set  absolutely,  using  an  inclinometer,  for  if  =  +90°  with 
the  quadrant  arm  horizontal. 

Lovell^  gives  a  detailed  description  of  the  calibration  scheme  used  for  the  probe, 
including  the  relationships  used  for  transforming  the  9,$  measurements  to  corrected 
flow  angles  in  pitch  (w°)  and  yaw  (v°)  -  and  also  gives  an  assessment  of  the  angular 
resolution  of  the  probe  and  the  precision  of  measurements.  In  essence  the  probe  calibra¬ 
tion  requires  the  following  quantities  as  functions  of  the  local  dynamic  pressure  q^ 

measured  by  the  probe: 

(i)  qj/q^  »  where  qj  is  the  true  dynamic  pressure  (Fig  9). 

(ii)  The  error  (H  -H.)/q,  where  H  is  the  total  head  measured  by  the  probe’s 

m  1  I  m 

central  tube  and  the  true  value  (Fig  10). 

The  calibrations  are  based  on  the  empty  tunnel  conditions  (ie  qt  =  q^,  H(  =  HQ,  in  this 
case),  the  'true'  values  being  established  by  measurements  with  a  pitot-static  tube.  The 
q^  date  shown  has  been  normalised  using  q^^  *  t*le  ^m  va*ue  at  the  &0  m/ s  test  speed. 
The  quantities  q^  and  have  been  measured  with  the  probe  nulled,  q^  is  define!  by 


qm  =  Hm  -  HP,  +  P2  +  P3  +  P4> 


where  Pj  ...  p^  are  the  readings  of  the  probes'  four  side  tubes.  Since  pressure  trans¬ 
ducers  of  adequate  quality  were  available  for  only  four  of  the  five  tubes,  the  calibration 
in  the  prevent  investigation  has  been  based  on 

=  Hm  "  ‘(PI  +  P2  +  2p3}  (3) 


where  each  of  the  four  pressures  is  measured  relative  to  the  tunnel  plenum  wall  static 
pressure  p^  .  The  probe  error  signals,  Apq  and  Ap^  ,  are  employed  in  the  probe  null¬ 
ing  procedure;  the  angles  9  and  <j>  are  adjusted  such  that  at  a  'balance'  condition 


ApQ  =  P1  "  p2  =  0 

(4) 

ApR  =  p3  "  p4  =  0 

(5) 

A  phase-advance  system  is  employed  to  increase  the  effective  frequency  response  of  the 
system,  by  compensating  for  the  viscous  damping  in  the  long  lines  between  the  probe  and 
the  transducers,  which  are  located  in  thermally  insulated  packages  outside  the  tunnel. 
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The  pressures  Ap^  and  Ap^  are  measured  independently  of  Hm,pj,p9,p^  ,  using  trans¬ 
ducers  functioning  as  part  of  the  probe  nulling  control  system.  (In  the  former  case 
transducers  with  minimal  zero  drift  are  required,  the  voltage  sensitivity  6 h / 6 p  being 
less  important;  in  the  latter  case  6E/6p  constancy  is  of  vital  importance  and  zero 
drift  effects  can  be  accommodated  by  keeping  measurement  runs  fairly  short.) 

The  calibration  scheme  also  requires: 

(i)  tp  as  defined  in  equation  (2)  above 

(ii)  the  pitch  of  the  traverse  gear  roll  axis  (a^) 

(iii)  the  yaw  of  the  traverse  gear  yaw  axis  (f^)  . 

Both  the  latter  are  referred  to  the  horizontal  centreline  of  the  tunnel  working  section. 
They  are  obtained  by  comparison  of  five-hole  probe  measurements  with  Conrad  tube  measure¬ 
ments  at  points  in  the  model  flowfield  where  substantially  pure  pitch  and  pure  yaw  con¬ 
ditions  exist. 

In  the  central  portion  of  each  vortex  core  the  normal  quadrant  angle  range 
(-30°  ^9  <  +30°)  (probe  located  at  the  quadrant  arm  centre)  has  been  found  insufficient 


to  achieve  probe  nulling  owing  to  extreme  values  of  the  flow  angles  (V°  —  W° 

max  max 


30.5°). 


Therefore  the  probe  was  remounted  at  one  extremity  of  the  60  circular  arc  quadrant  arm, 
and  the  range  of  quadrant  angle  0°  <  6  <  60°  then  readily  accommodated  the  outstanding 
measurements.  Some  earlier  measurements  (including  some  at  the  check  datum  point)  were 
repeated,  to  allow  the  re-evaluation  of  t  by  a  process  of  comparison. 

4. 5  Make  traversing  -  experiment  control,  data  logging,  data  reduction 

Ref  6  may  be  consulted  for  details  of  the  very  extensive  computer  software  available 
for  various  modes  of  control  of  wake  traverse  experiments,  for  logging  of  data,  and  for 
the  off-line  reduction  of  the  stored  data  (in  disk  file  form)  on  completion  of  measure¬ 
ment  runs. 

It  is  however  appropriate  to  note  the  following  points  here: 

(i)  Having  achieved  a  null  condition  the  analogue  voltages  representing  the 
measured  pressures,  the  current  orientation  (© , 4> )  of  the  probe  and  its  tip  position 
(x,y,z)  are  logged  by  the  computer  via  its  ADC  unit.  The  unprocessed  data  is  written  to 
a  disk  file  and  concurrently  output  as  a  single  line  of  line-printer  data.  At  the  end 
of  each  prescribed  sequence  of  measurements  the  corresponding  series  of  0,4>,pm  -  H 
datapoints  can  (if  required)  be  displayed  as  a  function  of  the  spatial  co-ordinate  being 
varied  (y  or  z  in  our  case)  on  an  analogue  graph-plotter. 

(ii)  Control  modes  exist  whereby  the  traverse  gear  will  make  measurements  at  a 
series  of  points  defining  lines  and  grids  on  the  traverse  plane  -  either  under  full 
computer  control  or  under  interactive  control  where  the  operator  guides  the  nulling 
process . 

(iii)  The  quality  of  the  raw  data  has  been  assessed  by  monitoring  the  variation  of 
the  voltages  corresponding  to  8 , d> , pm  -  H  during  each  run,  and  by  checking  the  consis¬ 
tency  of  measurements  at  the  check  datum  point  before  and  after  each  run. 
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(iv)  The  0,4i  nulling  operation  is  controlled  by  an  open-loop  control  system. 

The  current  6  and  $  values  at  any  instant  are  recorded  as  the  cumulative  sums  of  the 
starting  values  and  all  control  signals  (A0,A<J>)  subsequently  sent  to  the  quadrant  and 
roll  drive  motors.  These  cumulative  sums,  of  signal  pulses  performing  the  dual  function 
of  control  and  power  supply,  also  control  the  analogue  ramp  voltage  generator  which  is 
monitored  by  the  computer.  Thus  if  a  fault  occurs  in  the  0  or  counter  or  ramp 

generator  systems,  or  if  the  drive  motors  fail  to  respond  to  the  signal  pulses:  the 
probe  nulls  correctly;  the  probe  pressures  are  measured  correctly;  and  the  6  and/or  $ 
values  logged  at  the  null  condition  are  in  error.  A  recurrent  source  of  faults  is  the 
system  of  transmission  of  the  0  and  $  nulling  dr ive/ control  signals  through  a  low 
voltage  brush  system.  The  open-loop  nature  of  the  control  system  is  a  major  weakness  of 
the  current  wake-traverse  system;  the  addition  of  a  system  to  measure  0  and  1 j>  directly 
would  obviate  the  need  to  repeat  many  measurements.  Errors  in  0  and  if  are  frequently 
relatively  small,  but  in  the  drag  measurement  context  it  is  necessary  to  repeat  even 
marginally  suspect  measurements.  In  the  present  test  some  10-15%  of  excess  data  had  to 
be  recorded,  to  produce  the  5100  useful  points  free  of  0  and  <j>  errors. 


Wake- traversing 


analysis  of  processed  data 


The  measurements  on  the  traverse  plane  have  been  used  to  calculate  the  spatial 
distributions  of  the  integrands  in  equations  (2)  and  (3)  in  Fig  6b  -  ie  the  distributions 
of  the  local  drag  coefficient  contributions  per  mm  (Cj^  and  Cp j ^ ^  “  on  that  plane.  In 
the  case  of  the  wing  wake  (Fig  6a,  analysis  regions  1-5)  the  overall  CDj  and  cDjj 
contributions  are  each  calculated  in  two  stages,  the  first  being  a  process  of  vertical 
stripwise  integration  yielding  the  spanwise  variation  of  Cpj(y)  and  Cdjj(v),  (Figs  11 
and  12)  the  second  being  a  process  of  spanwise  integration.  In  the  case  of  the  strut 
wakes  it  is  convenient  to  define  Cpj.(z)  (Fig  13)  and  CDjj(z)  ,  and  the  order  of  integ¬ 
ration  is  reversed. 


All  these  integration  procedures  are  incorporated  in  the  computer  software  avail¬ 
able,  trapezium  rule  techniques  being  employed.  However,  as  with  the  spanwise  integration 

+  1 

of  lift  (section  3.3),  the  definitive  value  J~  Cp ^ (y) d (2y/b)  =  Cp^TyTb  =  Cp^  =  0.02287 


for  the  wing  wake  has  been  derived  by  graphical  integration  of  large  scale  plots.  The 
1.7%  disparity  from  the  computer  integration  is  largely  attributable  to  the  omission  of 
a  small  number  of  anomalous  Cp^  datapoints  from  the  analysis.  The  spanwise  mean  value 
of  drag  coefficient  contribution  per  unit  run  for  the  wing  wake  ^Cp^  ^  has  been  used  to 
normalise  the  profile  drag  data  presented  in  Figs  11  and  13. 


The  spanwise  integration  of  the  vortex  drag  data  has  been  performed  using  a  computer 
procedure.  Fig  12  shows  the  spanwise  variation  of  Cpjj(y)  ;  again  the  data  have  been 
normalised  using  a  spanwise  mean  value  ^Cp^j(y)b  =  0.00937^  appropriate  to  the  wing 
wake  (Fig  6a,  analysis  regions  1-5). 


Maskell's  theory,  in  the  form  presented  in  Ref  1,  assumes  that  on  the  upstream 
reference  plane  the  centreline  reference  total  head  H()  is  valid  over  the  whole  plane. 
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Examination  oi  the  present  data  shows  that  slightly  negative  total  head  defects 
(AH/q()  =  (H()  -  H )  /  q  ( ( )  have  been  recorded  outside  the  viscous  wake,  particularly  towards 
the  wing  tips.  Small  adjustments  have  been  made  in  the  AH/q(J  data  (before  calculation 
of  the  flow  velocity  components  from  the  V*,W°  ,  and  static  pressure  data)  to  compensate 
for  this.  A  representative  mean  correction  has  been  applied  in  each  of  eleven  spanwise 
wing  wake  segments  and  in  the  strut  wake  regions,  the  correction  being  at  most 
AH/q^  =  0.0017  .  These  adjustments  correct  for  the  small  spanwise  variation  of  free- 
stream  total  head  from  the  centreline  reference  value  Hf)  -  and  in  fact  correspond  to 
the  spanwise  variation  of  H  recorded  in  the  empty  tunnel  (Fig  14,  z  =  0)  .  These 
AH/q^  corrections  amount  to  17-  of  the  final  value.  No  similar  correction  is 

appropriate  to  the  Cp^  and  vorticity  calculations. 

5  EXPERIMENTAL  RESULTS  -  GENERAL  COMMENTS 

5.  1  Forces  and  moments  data 

Figs  3,  4a  and  5  show  the  basic  lift,  drag  and  pitching  moment  characteristics  for 
both  transit ion-f ixed  and  transit ion-f ree  conditions.  The  mean  and  values 

recorded  during  wake-traversing  (transition  fixed)  are  in  good  agreement  with  the  approp¬ 
riate  curves.  All  this  data  derives  from  mechanical  balance  measurements,  and  it  has 

9  ... 

been  corrected  using  conventional  techniques.  The  support  rig  correction  is  restricted 

to  the  deduction  of  the  model-off  tares  measurements  (section  4.3).  Correction  for 

tunnel  constraint  effects  on  the  drag  polar  has  been  made  using  the  Glauert  expression 

ACDi  =  (6) 

which  is  considered  adequate  for  the  present  moderate  lift  application. 

Fig  4b  shows  the  low  lift  segment  of  the  drag  polar,  which  has  a  minimum  drag  point 
(0^^  Cpc)  at  =  f'jj  =  0.0203  .  In  a  restricted  range  of  lift  values 

(0.3  <  CL  <  0.6)  encompassing  the  traverse  condition  (Cl  =  0.491),  a  lift-dependent 

C  v  c 

drag  factor  value  k  =  1.42  holds  for  the  above  C  ,  where  k  is  defined  by  assumption 

1 4  U 

of  a  parabolic  lift-dependence  of  the  total  drag 
5 . 2  General  characteristics  of  the  flowfield 

The  principal  features  of  the  wake  are  shown  by  the  contours  of  constant  total  head 
defect  (AH/q(1)  measured  in  the  traverse  plane,  Figs  15a  and  16  respectively  showing  the 
wake  behind  the  fuselage  and  behind  the  port  wing.  Fig  17  shows  the  locus  of  maximum 
AH/qfl  in  the  viscous  wake  from  the  port  wing  and  the  spanwise  variation  of  maximum  total 

head  defect  AH/q()  .  At  all  points  on  the  locus  AH/q^  is  very  much  less  than  the 

values  observed  at  the  centres  of  the  tip  vortex  cores  (AH/q^  =*  1.4)  and  considerably 
less  than  the  maximum  value  observed  in  the  fuselage  wake  (AH/q^  —0.7). 

The  points  A,  A'  and  B,  B'  marked  on  Fig  15a  represent  the  estimated  locations  of 
the  vorticity  maxima  (Fig  15b)  produced  by  the  two  pairs  of  vortices  assumed  to  originate 
in  the  wing- f use lage  junctions,  and  these  four  points  lie  in  a  region  of  moderately  high 
',H/q  .  The  contours  of  highest  AH/q^  are  closely  consistent  with  the  marked  locations 

(y,7.)  of  the  fuselage  tailcone  tip  ('T',  Fig  1)  and  the  pitch  strut  pivot.  The  large 
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area  of  significant  AH/qg  below  the  fuselage  wake  is  largely  attributable  to  an 
unfavourable  interaction  between  the  pitch  strut  and  the  fuselage.  Fig  15a  also  shows 
•■he  projections  of  the  fuselage  cross-sections  at  station  x  »  0  (the  plane  of  the 
trai l ing-edge  tip-rip  line)  and  x  =  -763  mm  (the  plane  of  the  wing  planform  apex)  onto 
the  traverse  plane. 

Figs  15a,  16  and  17  also  show  that  the  wing  wake  is  significantly  thinner  in  the 
vicinity  of  the  wing  root  and  tip  than  in  the  mid  semispan  region,  and  that  this  thinning 
is  accompanied  by  a  reduction  in  maximum  AH/qg  .  The  maximum  AH/q0  falls  towards 
zero  as  one  follows  the  residual  viscous  wake  round  the  circumference  of  the  tip  vortex 
core . 

Figs  15b  and  18  show  the  corresponding  contours  of  constant  local  streamwise 
vorticity  £  (in  nondimensional  form)  in  the  traverse  plane  aft  of  the  fuselage  and  the 
port  wing.  Fig  19  shows  the  z  value  at  which  maximum  (negative)  vorticity  was  recorded 
at  each  spanwise  station  considered,  together  with  the  spanwise  variation  of  this  maximum. 
Disregarding  the  wing  root  region,  where  the  £  values  were  so  small  as  to  make  it 
difficult  to  determine  maxima,  the  z  positions  of  the  £  maxima  only  deviate  signific¬ 
antly  from  the  locus  of  maximum  AH/qg  (Fig  17)  in  the  strut/wing  intersection  region 

(2y/b  -  0.45),  The  maximum  £  value  increases  spanwise  from  root  to  tip,  until  the 
viscous  wake  becomes  indistinguishable  from  the  vortex  core  in  £  terms.  A  localised 
reduction  in  £  occurs  in  the  strut  top  region  just  above  the  point  where  vorticity  is 
generated  that  is  of  opposite  sign  to  that  in  the  adjacent  tip  vortex. 

Four  pairs  of  vortices  are  discernible  in  the  centre  section  wake  (Fig  15b),  all 

regions  of  significant  vorticity  being  enveloped  by  the  AH/qQ  =  0.025  total  head  defect 

contour.  Regions  A  and  A1,  and  B  and  B',  are  clearly  associated  with  vortices  respec¬ 
tively  generated  in  the  upper  and  lower  wing  roots;  regions  C  and  C'  appear  to  be  associ¬ 
ated  with  a  vortex  pair  springing  from  the  forward  fuselage  roof,  though  their  precise 
origin  is  unclear;  finally,  regions  D  and  D1  are  associated  with  a  vortex  pair  generated 
by  the  sharp  inclined  upper  edges  of  the  top  of  the  pitch  strut.  The  contours  of  non- 
dimensional  vorticity  of  magnitude  |£c/Ug[  =  0.05  represent  the  limit  of  vorticity 
resolution  with  the  wake  traverse  system  in  this  particular  flowfield. 

Fig  18  shows  the  contours  of  vorticity  for  the  port  wing  wake  and  their  relation¬ 
ship  to  the  contour  AH/qg  =  0  which  defines  the  estimated  'edge'  of  the  viscous  wake. 

The  strut/wing  intersection  secondary  flow,  containing  vorticity  of  opposite  sense  to 
that  in  the  tip  vortex,  is  a  prominent  feature.  All  essential  features  shown  in  the 
figure  are  mirrored  in  the  starboard  wing  wake. 

The  circulation  y  has  been  evaluated  on  rectangular  contours  in  the  traverse 
plane  enclosing  the  centres  of  the  eight  fuselage-region  vortices,  the  strut  secondary 
flow  vortices  and  the  tip  vortices  -  and  the  results  have  been  expressed  in  a  form 
analagous  to  that  of  lift  coefficient  increments,  as  shown  in  the  following  table. 

Although  the  fuselage  vortices  are  individually  quite  strong  relative  to  the  tip  vortices, 
their  net  contribution  to  the  vortex  drag  of  the  configuration  (Fig  6b,  equation  (3))  is 
negligible,  a  result  accentuated  by  the  weighting  effect  of  the  small  V'(y,z)  values  in 
the  centre  section. 
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Nondimensional  circulation  (l'Y/Uf)c)  of  vortices 


Vortex 

Fort 

Starboard 

Tip  vortex 

-0. 280 

+0. 271 

Strut/wing 

intersection  vortex 

+0.014 

-0.  01  6 

A’  ,  A 

+0.065 

-0.04') 

Fuselage 

B',  B 

-0.01 1 

+0 . 0 1 A 

vort ices 

c' ,  C 

-0.007 

+0.007 

o'  ,  D 

-0.026 

+0.013 

Figs  20  and  21  respectively  show  plots  of  the  contours  of  constant  sidewash  angle 
and  constant  downwash  angle  defined  by  the  wake-traverse  measurements.  In  both  cases  a 
pronounced  increase  in  the  magnitude  of  local  flow  angles  occurs  in  the  vicinity  of  the 
strut/wing  intersection.  A  noteworthy  feature  of  the  latter  figure  is  the  strong  upwash 
region  that  occurs  on  the  centreline  (y  =  0)  behind  the  fuselage.  This  is  consistent  with 
the  existence  ot  the  strong,  upper  wing-root  junction  vortices  A  and  A1  inferred  from  the 
vorticity  contours  ('Fig  lib),  such  chat  fluid  tends  to  flow  from  the  upper  wing  surfaces 
up  the  fuselage  sides,  thereby  thinning  the  wing  boundary  layer  and  the  consequent  wake 
in  proximity  to  the  fuselage. 

Figs  22  and  21  show  the  downwash  angle  and  total  head  defect  (AH/q^)  data  produced 
by  limited  measurements  close  to  the  fuselage  side,  at  stations  ’A’  and  ' B ’  respectively 
(Fig  1).  These  measurements  were  made  in  an  attempt  to  trace  the  path  of  the  upper  wing- 
root  junction  vortex  A  (Fig  15b)  along  the  fuselage  side  -  in  the  expectation  that  its 
presence  would  produce  a  significant  distortion  of  the  fuselage  boundary  layer  and  locally 
increased  values  of  AH/ qQ  .  The  figures  show  this  expectation  to  have  been  unjustified, 
the  measurements  extending  a  good  way  into  the  fuselage  boundary  layer  with  o  fine  mesh 
(Ay  =  Az  =  4  mm).  The  downwash  data  in  Fig  22  gives  some  ambiguous  indication  of  a  fuse- 
Lage  side  vortex  at  station  'A'  -  at  station  ' B'  (Fig  21)  it  appears  that  the  vortex  has 
migrated  so  far  downwards  that  it  is  below  the  z  range  in  which  downwash  measurements 
were  made.  In  such  close  proximity  to  the  fuselage  it  was  not  possible  to  null  the  probe 
fully  {is  to  roll  the  quadrant  arm  from  the  vertical,  i,  =  &  ,  position)  and  to  determine 
the  sidewash  field.  Therefore  the  w°  and  AH/q^  calculations  depend  on  the  assumption 
of  small  sidewash  angles,  and  it  has  not  been  possible  to  determine  the  vorticity  field 
close  to  the  fuselage.  On  the  basis  of  this  limited  evidence  it  can  tentatively  be  infer¬ 
red  that  the  wing-root  junction  vortices  A  and  A'  remain  highly  concentrated  along  the 
parallel  part  of  the  fuselage,  and  migrate  downwards  in  the  downwash  field  generated  by 
the  wing.  Over  the  converging  rear  fuselage,  where  some  degree  of  adverse  pressure 
gradient  exists,  it  is  probable  that  the  vortices  A  and  A'  diffuse  rapidly  and  migrate 
rapidly  upwards  to  their  observed  positions  in  the  principal  wake-traverse  plane 
(Fig  1  )  ad b ) . 


Great  care  was  taken  to  set  the  model  up  symmetrically  with  respect  to  the  tunnel 
and  traverse  gear  centreline  plane  (y  =  0).  In  the  mean  the  net  crossflow  generated  by 
the  model  was  small;  an  analysis  of  720  pairs  of  points  symmetrically  disposed  to  port 
and  starboard  in  the  wake  (for  AH/q^  ^  0.05)  yielded  a  mean  yaw  of  only  0.06°.  It  has 
also  been  found  that  the  effective  aerodynamic  centreline  plane  of  the  model  closely 
agrees  with  the  geometric  centreline  plane.  Examination  of  the  locations  of  peak  local 
flow  angles,  peak  local  vortex  drag  contribution  (Cp^)  and  peak  local  profile  drag 
contribution  (Cp^ )  consistently  locate  the  port  and  starboard  vortex  core  centres  at 
spanwise  positions  2y/b  «=  -0.953,  +0.959  respectively.  A  similar  procedure  fixes  the 
heights  of  the  vortex  core  centres  as  z  =  20,  24  mm  for  the  port  and  starboard  vortices, 
respectively.  This  4  mm  height  discrepancy  (0.2%  of  the  span)  may  well  be  attributable 
to  small  asymmetries  in  model  dimensions  in  the  tip  regions. 

Figs  24  to  27  allow  further  assessment  of  the  extent  of  port-to-starboard  symmetry 
in  the  wake,  for  two  z  heights}  the  z  =  0  line  approximates  to  the  upper  edge  of  the 
wake  and  the  z  =  -80  mm  line  runs  some  little  way  below  it.  Figs  24  and  25  show  the 
spanwise  variation  of  downwash  angle  at  the  two  heights  in  the  traverse  plane,  while 
Figs  26  and  27  show  the  corresponding  sidewash  angle  data.  A  high  degree  of  port- 
starboard  symmetry  exists  over  the  bulk  of  the  span,  the  apparent  asymmetry  in  the  tip 
regions  being  attributed  to  the  above  mentioned  small  asymmetries  in  vortex  core  location 
relative  to  the  traverse  gear  axes. 

Figs  28a-c  show  the  variation  of  the  three  velocity  components  (U,v,w)  along 
vertical  and  horizontal  lines  through  the  port  tip  vortex  core  in  the  traverse  plane. 

For  the  horizontal  cuts  only  a  partial  set  of  experimental  data  is  available  for  the 
precise  vortex  centre  height  (z  =  +20  mm)  and  a  set  of  data  for  z  =  +24  mm  is  addition¬ 
ally  shown.  Fig  28d  shows  the  corresponding  variations  of  total  pressure  and  static 
pressure  defect  vertically  through  the  vortex  centre.  Fig  28e  shows  the  corresponding 
vertical  variation  of  local  streamwise  vorticity  (5);  in  this  case  some  data  from  the 
starboard  vortex  has  been  shown  for  comparison.  For  both  the  port  and  starboard  cuts 
the  £(z)  data  for  the  vortex  core  has  been  evaluated  using  mesh  sizes  of  both  4  x  4  mm 
and  8*8  mm.  The  coarse  mesh  data  generally  defines  satisfactory  local  £  values,  but 
the  definition  of  the  £(z)  curve  is  insufficient  to  define  the  maximum  at  the  vortex 
centre  adequately.  (An  analagous  exercise  on  the  spanwise  variation  of  the  local  vortex 
drag  contribution  per  unit  run  (Cp^iy))  showed  that  values  calculated  from  coarse  mesh 
(8  mm)  data  were  accurate  enough  in  themselves,  but  that  4  mm  data  were  required  to 
establish  the  maximum  satisfactorily.) 

Fig  29  shows  the  variation  of  downwash  angle  along  lines  of  constant  y  in  the 
traverse  plane,  through  the  port  and  starboard  vortex  centres,  and  Fig  30  shows  the 
corresponding  variation  of  sidewash  angle.  The  port-starboard  disparities  in  the  flow 
angles  at  stations  2y/b =  -0.953,  +0.959  may  result  from  these  cuts  not  passing  exactly 
through  the  vortex  core  centres.  Downwash  data  have  been  plotted  on  Fig  29  for  y  cuts 
0.6%  semispan  outboard  of  and  0.4%  inboard  of  the  cut  most  nearly  passing  through  the 
port  vortex  core  centre,  to  demonstrate  the  large  spanwise  variation  of  downwash  profile 
in  the  core  region. 


Fig  31  shows  the  contours  of  constant  total  head  defect  (AH/q())  established  in  a 
survey  ot  the  port  strut  guard  leakage  plume  in  the  principal  traverse  plane;  a  similar 
picture  was  obtained  on  the  starboard  side.  At  its  minimum  z  p-sition  (-800  mm  relative 
to  the  current  tip-tip  line)  the  traverse  gear  indicated  the  upper  edge  of  a  rather 
smaller  plume  from  the  pitch  strut  guard.  These  plumes  are  produced  by  the  spillage  of 
low  energy  tluid  over  the  top  of  the  strut  guards;  this  fluid  originates  outside  the 
tunnel  and  is  drawn  up  the  gaps  between  the  struts  and  their  windshields  (or  guards)  by 
the  sub-atmospheric  pressure  level  in  the  working  section.  Clearly  no  significant 
portion  of  the  momentum  defect  in  these  plumes  is  attributable  to  drag  on  the  exposed 
portions  of  the  struts. 

6  EXPERIMENTAL  RESULTS  -  WAKE- TRAVERSE  DRAG 
6 . 1  General  comments 

Fig  11 a&b  show  the  spanwise  variation  of  local  contribution  from  the  wing  and 

fuselage  wakes,  while  Fig  12a-c  show  the  corresponding  Cp^  data.  Fig  1 3a&b  show  the 
variation  of  Cd^  contribution  along  the  main  struts  and  the  pitch  strut  respectively, 
there  being  no  perceptible  Cp,^  contribution  from  the  strut  wakes.  All  this  data  is 
normalised  with  respect  to  the  spanwise  mean  values  of  local  Cp^  or  cDj j  contribution, 
these  means  being  calculated  on  the  basis  of  exclusion  of  strut  drag  contributions  and 
of  inclusion  of  the  estimated  strut-on-model  interference  effects. 

The  wake  of  the  model  and  its  support  rig  have  been  divided  into  arbitrary  analysis 
regions  (Fig  6a,  Table  A).  Inevitably  the  drag  increments  measured  in  the  fuselage  and 
wing  wake  regions  will  include  contr ibut ions  from  short  lengths  of  the  three  struts,  and 
from  the  viscous  interference  effects  of  the  struts  on  the  model.  The  strut  wake  analysis 
regions  have  been  specified  so  as  to  exclude  all  fuselage  and  wing  wake  drag  contributions, 
together  with  those  from  the  main  strut/wing  junction  secondary  flows.  They  include  the 
bulk  of  the  model-on-strut  interference  effects  and  all  contributions  from  the  readily 
identifiable  portions  of  the  upper  strut  wakes.  The  small  drag  contributions  from  the 
three  unrecessed  cleat  fittings  have  been  assumed  negligible  -  the  effects  of  the  pitch- 
strut  cleat  being  included  in  the  fuselage  drag,  those  from  the  main  strut  cleats  being 
included  in  the  strut-on-model  interference  terms. 

Figs  II  and  12  show  that  a  high  degree  of  port-to-starboard  symmetry  exists  in  the 
measured  spanwise  distributions  of  local  Cp^  and  Cp^  contributions.  This  symmetry 
is  reflected  in  the  integrated  drag  contributions  for  the  various  wake  analysis  regions 
shown  in  Table  4a&b.  Section  6.3  details  the  means  by  which  the  model-alone,  interference 
of  support  rig  on  model,  and  support  rig  (net  of  interference  of  model)  drag  increments 
have  been  isolated.  Sections  6. A  and  6.5  detai1  the  procedures  used  to  apply  tunnel 
constraint  corrections  and  to  establish  the  induced  drag  from  the  corrected  ^Djj  data. 

Examination  of  Figs  lib  and  12b  shows  that  just  outboard  of  the  main  strut  locations 

a  large  perturbation  from  the  general  trend  occurs  in  both  the  Cp^.  and  cDji  data. 

The  double-peak  feature  on  the  Cp^  distribution  has  been  attributed  partly  to  a  contrib¬ 
ution  from  the  short  length  of  main  strut  unavoidably  included  in  the  wing  wake  analysis 
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region  (inboard  peak),  and  partly  to  the  secondary  flow  generated  at  the  strut/wing 
intersection  (outboard  peak)  (see  also  Fig  18). 

The  spatial  distributions  of  profile  drag  (Cjjj(y,z))  and  vortex  drag  (C£jjj(y,z)) 
contribution  on  the  traverse  plane  may  be  inferred  from  the  contour  plots  of  total  head 
defect  (Fig  16)  and  streamwise  vorticity  (Fig  18)  -  together  with  Figs  11  and  12  -  and  it 
has  not  been  considered  necessary  to  include  contour  plots  of  CjjT  and  in  this 

Report. 


Comparison  of  Fig  12a  -  showing  the  spanwise  variation  of  in  the  fuselage 

and  wing-root  region  -  with  Fig  15b  -  which  shows  a  contour  plot  of  streamwise  vorticity 
(closely  similar  to  Cp^iy.z))  in  the  same  region  -  illustrates  the  relatively  limited 
utility  of  the  former  presentation  of  data  for  flowfield  diagnostic  purposes.  A  net 
Cqjj  of  virtually  zero  at  prescribed  y  can  in  fact  be  the  resultant  effect  of  a  series 
of  vorticity  concentrations  centred  at  various  z  positions  at  that  y  . 

6.2  Comparison  with  mechanical  balance  data 

The  most  obvious  criterion  of  the  accuracy  of  the  wake-traverse  technique  is  that 
used  by  Maskell1  -  namely  a  comparison  with  mechanical  balance  measurements.  Clearly  such 
a  comparison  should  be  made  on  the  basis  of  the  respective  uncorrected  lift  and  drag 
figures  for  the  model  together  with  that  portion  of  its  support  rig  exposed  to  the  tunnel 
stream.  This  obviates  the  two  subjective  judgments,  possibly  inconsistent  ones,  required 
to  extract  the  respective  support  rig  and  interference  effects.  The  essential  remaining 
assumption  is  that  both  the  traverse  system  and  the  balance  system  are  sensing  the  same 
effect.  The  validity  of  this  is  by  no  means  so  clear  for  the  present  configuration  as 
for  Maskell 's  wire  support  rig. 


The  traverse-derived  drag  coefficient  is  Cpt  =  0.03875  ,  which  includes  the  three 
auxiliary  terms  (ACp  ,  ACdtt  ,  ACqtt  ;  section  3.2).  In  those  regions  of  the  strut  wakes 
where  strut  effects  are  indistinguishable  from  those  of  the  low  momentum  strut  guard 
leakage  plumes  (Fig  31),  the  strut  drag  contributions  have  been  estimated  by  extrapola¬ 
tion  of  the  mean  Cp^(z)  levels  in  the  central  portions  of  the  struts  to  the  planes  of 
the  strut  guard  tops.  This  assumes  that  the  leakage  plumes  have  no  material  effect  on 
the  Cpj(z)  distributions  that  would  be  discernible  in  their  absence.  On  this  basis  of 
comparison  the  traverse-derived  drag  coefficient  exceeds  the  apparently  equivalent 
balance-derived  value  (0.03866)  by  0.2%. 


As  a  further  check  on  the  general  accuracy  of  the  traverse  technique  it  was  decided 
to  compare  traverse-derived  and  balance-derived  drag  data  for  the  two  main  struts  in  a 
model-off  condition.  (It  was  considered  justifiable  to  double  the  traverse-derived 
measurement  of  the  drag  of  a  single  strut  (the  starboard  one,  in  fact)  and  to  compare 
this  with  the  balance-derived  figure  for  an  identical  pair  of  struts.)  The  traverse 
measurements  were  made  quite  close  to  the  strut  (x  =  -96  mm),  securing  advantages  of 
measurement  accuracy  in  a  compact  wake  and  leakage  plume,  and  the  results  are  shown  in 
Fig  13a(i)&(li)  (together  with  model-on  data,  to  illustrate  the  strong  negative  model-on- 
struts  interference).  The  traverse-derived  strut  drag  was  calculated  by  graphical  integ¬ 
ration  as  ACpj  =  0.00217  ,  again  using  the  Cpj(z)  extrapolation  procedure  in  the  region 
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of  the  leakage  plume.  The  corresponding  mechanical  balance  figure  was  AC[)j  =  0.00284. 

(In  an  empty  tunnel  the  struts  are  non-lifting  bodies  and  no  appreciable  CDjj  effects 
are  generated.) 

The  results  of  this  comparison  -  giving  a  ACq^  =  0.00067  discrepancy  for  the  single 
main  strut  -  are  shown  in  Table  5a;  the  traverse-derived  drag  increments  relate  to  the 
same  (y,z)  analysis  region  bounds  as  used  in  model-on  traverse  measurements.  This  result 
makes  it  apparent  that  the  two  systems  are  not  in  fact  sensing  the  same  effects.  It  has 
been  concluded  that  the  balance  registered  a  ACp  =  0.00067  contribution  attributable  to 
buoyancy  effects  on  the  shielded  portion  of  the  strut  within  the  strut  guard.  This  effect 
is  analagous  to  a  well-established  buoyancy  effect  on  lift  tares  measurements 1 ^  -  both 
effects  being  thought  to  be  produced  by  the  leakage  flow  into  the  working  section  through 
the  clearance  passages  between  struts  and  guards.  Obviously  a  momentum  defect  relating 
to  the  buoyancy  drag  increment  will  appear  in  the  wake-traverse  plane  but,  like  the 
momentum  defect  attributable  to  the  lower  unshielded  portions  of  the  struts,  will  be 
indistinguishable  from  the  main  leakage  plume  effects  (Fig  31).  It  is  not  considered 
possible  that  any  significant  portion  of  the  ACp  =  0.00067  disparity  is  attributable  to 
deficiencies  in  the  Cj^fz)  extrapolation  assumptions. 

It  should  be  noted  that  the  tares  measurements  (balance),  model-off  traverse 
measurements  and  model-on  traverse  measurements  were  made  at  different  times,  that  the 
relative  positions  of  the  struts  and  their  guards  were  likely  to  have  varied  (in  ignorance 
of  the  significance  of  this  factor)  and  that  the  ACp  =  0.00067  buoyancy  effect  estimate 
is  therefore  subject  to  a  measure  of  uncertainty.  It  was  not  practical  to  undertake  a 
similar  exercise  for  the  pitch  strut,  but  a  reasonable  estimate  for  the  buoyancy  drag 
increment  for  this  would  be  half  the  above  figure. 

In  the  light  of  the  above  results  it  is  clear  that  an  acceptable  basis  of  comparison 
between  the  traverse-derived  drag  of  the  model  and  its  support  rig  (which  involves  an 
extrapolation  to  eliminate  all  leakage  flow  related  effects)  and  the  balance-derived  drag, 
requires  the  latter  to  be  net  of  a  AC^  =  -O.OOI67  correction  for  the  leakage  buoyancy 
effects  for  the  three  struts.  On  this  basis  the  wake-traverse-derived  drag  (Cp  =  0.03875) 
exceeds  the  corresponding  balance-derived  figure  (Cp  =  0.03699)  by  4.8%. 

Table  2  shows  a  comparison  between  balance  and  traverse  derived  lift  and  drag 
figures  for  the  current  investigation,  and  for  two  others^’^  on  a  model  (BAe  146)  with 
74%  of  the  current  planform  reference  area.  The  above  process  of  deduction  of  leakage 
flow  effects  has  been  consistently  applied,  and  in  each  case  the  data  refers  to  the  model 
in  combination  with  its  support  rig.  The  balance-traverse  drag  disparities  are  worse  in 
percentage  terms  than  for  the  present  model  -  but  they  are  more  closely  comparable  when 
expressed  in  absolute  terms  and  scaled  on  a  common  reference  area  (that  of  the  current 
model),  as  follows: 
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Model 

Ct 

UC 

ACDabs 
( t  rav-bal ) 

AC 

( t  rav 

-bal ) 

(a) 

(b) 

Current 

0.49 

+0.00176 

+  0.2 

+  4.8 

BAe  146 

0.  30 

+0.00232 

+  2.5 

+  9.4 

BAe  146 

0.61 

+0.00208 

+  3.0 

+  6.5 

The  above  tabulated  data  suggest  that  one  has,  within  the  uncertainty  of  the 
measurements,  a  virtually  constant  absolute  disparity  between  the  balance-derived  and 
traverse-derived  results.  This  would  be  consistent  with  an  error  in  the  traverse 
technique  originating  in  a  support  rig  effect.  It  is  important  to  note  that  the  same 
support  rig  has  been  employed  in  all  three  cases  (though  possibly  with  some  variation  of 
strut-guard  clearance  passage  characteristics)  -  and  that  the  one  set  of  model-off  strut 
traverse  measurements  has  been  used,  together  with  each  of  the  three  model-off  strut 
tares  balance  measurements,  to  estimate  the  leakage  buoyancy  effect  in  each  case.  Any 
deficiency  in  model-off  strut  traverse  data  would  therefore  be  reflected  as  a  consistent 
absolute  ba lance/ traverse  disparity  in  each  of  the  three  sets  of  data.  Columns  (a)  and 
(b)  of  the  above  table  correspond  to  the  percentage  drag  disparities  between  balance  and 
traverse  data,  according  as  to  whether  the  leakage  buoyancy  effect  as  quantified  is 
respectively  neglected  or  allowed  for.  A  careful  check  of  the  model-off  strut  traverse 
data  suggests  that  the  latter  is  the  correct  approach,  and  that  there  is  some  consistent 
source  of  error  in  the  technique  as  currently  applied. 

6 . 3  Extraction  of  support  rig  effects 

Table  4  shows  the  breakdown  of  total  traverse-derived  drag  coefficient  (Cp^)  into 
CDl  and  contributions  relating  to  the  wing  and  fuselage  wake  (analysis  regions 

1  —  5 )  and  to  the  support  rig  (analysis  regions  6-8).  No  attempt  has  been  made  to  distri¬ 

bute  similarly  the  small  net  contribution  from  the  three  auxiliary  terms  (AC^  =  -0.00053). 
The  wing  and  fuselage  increments  include  the  estimated  strut-on-model  interference  effects 
of  the  two  main  struts  within  the  =  0.03224  total.  In  the  present  simple  experiment 
it  has  not  been  possible  to  quantify  the  strut-on-model  interference  contribution  of  the 
pitchstrut,  which  is  likely  to  be  very  small,  and  this  has  been  assumed  to  be  zero.  (The 

drag  attributable  to  the  pitchstrut  cleat  has  been  accounted  to  the  fuselage,  and  that  of 

the  main  strut  cleats  to  the  strut-on-model  interference  terms.)  The  support  rig  drag 
increment  (C^  =  0.00704)  comprises  the  model-off  strut  drag  (less  an  allowance  for  the 
effects  of  the  exposed  strut  tops)  and  the  model-on-strut  interference  effects.  It  is 
not  possible  directly  to  identify  the  drag  increments  attributable  to  the  short  lengths 
of  the  three  struts  projecting  into  the  model  wake  analysis  regions.  These  increments 
have  been  estimated  by  extrapolation  to  the  strut  tops  (Fig  1 3a(i)fib (i) ) . 
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The  CDl  profile  drag  interference  terms  in  Table  4a  have  been  estimated  by 

fairing-out  the  strut  region  perturbations  in  the  port  and  starboard  Cpj(y)  distribu¬ 
tions  (Fig  I  I  a) ,  so  establishing  the  sum  of  the  upper  strut  contributions  and  interference 
contributions.  Having  estimated  the  former  by  extrapolation  the  latter  may  be  deduced, 
and  in  the  present  case  the  two  effects  are  of  broadly  similar  magnitude.  The  inter¬ 
ference  term  quantifies  the  disturbance  to  the  wing  pressure  field  and  undersurface 
boundary  layer  produced  by  the  strut  (and  its  cleat  fitting)  and  the  secondary  flow 
generated  by  the  wing/strut  junction  (Fig  18),  so  increasing  . 

The  Cn . T  .  vortex  drag  interference  terms  in  Table  4(a)  have  been  estimated  by  a 
11  mt 

similar  curve-fairing  exercise  on  the  (y)  distributions  (Fig  12b).  Here  a  much 

more  subjective  judgment  is  required,  to  assess  the  mean  spanwise  trend  of  Cj)^  in  the 
absence  of  a  strut.  As  with  the  analysis  of  Co^(y)  one  is  obliged  to  assume  that 
observable  local  perturbations  represent  the  entire  net  interference  effect.  In  fact  the 
localised  modification  of  shed  vorticity  will  affect  the  whole  spanwise  load  distribution, 
and  there  will  be  second-order  interference  effects  on  the  distribution  in  addition 

to  the  first-order  effects  one  is  able  to  identify  and  correct. 

Table  5(a)  shows  the  process  of  estimation  of  the  leakage  flow  buoyancy  drag  incre¬ 
ment  for  the  main  struts,  discussed  in  section  6.2.  Table  5(b)  shows  the  build-up  of  the 
support  rig  drag  contributions  (Cp  =  0.00704  in  total)  entered  in  Table  4(b).  It  is 
possible  to  estimate  the  model-on-strut  interference  drags  in  Table  5(b)  by  comparison  of 
the  model-on  drag  contributions  (traverse  data)  with  the  model-off  cares  data  (mechanical 
balance)  -  having  allowed  for  the  effects  of  the  leakage  flow  (the  pitch  strut  value 
being  taken  as  half  that  for  the  main  struts)  and  the  exposed  strut  tops.  In  the  case 
of  the  main  struts  this  interference  is  negative,  as  is  qualitatively  apparent  from 
Fig  1  3a( i) . 

In  the  case  of  the  pitch  strut  one  estimates  a  positive  interference  drag  increment 
Cp  =  0.00098  .  The  Cp  =  0.00325  total  pitchstrut  drag  includes  two  contributions, 
denoted  'A'  (0.00050,  the  drag  extrapolation  for  the  strut  top)  and  ’B’  (0.00088)  on 
Fig  13b(i)  -  the  latter  being  readily  identifiable  as  an  adverse  interaction  between  the 
model  and  the  pitch  strut  (Fig  15a),  rather  than  a  fuselage  flowfield  effect.  In  fact 
a  pair  of  vortices  (Fig  15b;  denoted  D  -  D ’ )  originates  in  this  region,  probably  generated 
by  separations  from  the  sharp  edges  of  the  oblique,  square-cut,  aft-sloping  top  end  of 
the  aerofoil  section  pitchstrut  (whose  oblique  end  plane  runs  into  the  t rai 1 ing-edge  line 
at  z  -  -70  mm  ).  The  extent  of  consistency  between  the  inferred  interference  increment 
(Cp  =  0.00098)  and  the  drag  increment  'B'  (Cp  =  0.00088)  suggests  firstly  that  the  above 
separation  effects  are  insignificant  in  the  model-off  case,  and  secondly  that  the  Cp 
increments  assumed  for  the  exposed  strut  top  effect  (0.000)2)  and  the  leakage  effect 
(0.00033)  were  not  seriously  unrealistic. 

6.4  Tunnel  corrections 

Table  3  shows  a  comparison  between  the  total  corrected  mechanical  balance  drag 

(Cn  J  and  the  equivalent  traverse  figure  (Cy.  ).  The  corrections  for  tunnel  con- 
cbal  C?rP 

straint  and  stream  pitch  ('E'),  solid  blockage  (  F  ),  empty  tunnel  streamwise  buoyancy 
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('G'),  and  the  omission  of  a  correction  for  support  rig-on-model  effects  are  common  to 
both  sets  of  data.  A  disparity  of  treatment  of  data  inevitably  arises  for  the  model-on- 
rig  effects,  and  for  this  reason  the  comparison  is  not  as  meaningful  or  satisfactory  as 
that  in  Table  2  (see  section  6.2).  The  deficiencies  of  the  Table  3  comparison  show  that 
a  satisfactory  evaluation  of  traverse-derived  net  model  drag  data  in  terms  of  balance- 
derived  data  is  possible  only  if  the  latter  includes  comprehensive  mechanical  balance 
tares  and  interference  measurements. 

b . 5  Induced  drag 

This  has  been  calculated  from  the  traverse  data  as  shown  in  the  right-hand  half  of 
Table  3.  The  measured  CDlI  term,  taken  literally  as  'vortex  drag'  for  the  present 
investigation,  is  0.00937  in  coefficient  form,  or  0.00889  net  of  the  two  auxiliary  terms 
(Table  4(c)).  It  has  been  argued  (section  3.1)  that  the  CDj.  contribution  measured  in 
the  vortex  cores  (.'.Cp  =  0.00044)  is  correctly  accountable  to  the  Cp  term,  and  this 

transfer  has  been  made.  A  correction  has  also  been  applied  for  stream  pitch  (-0.00043), 
and  the  Glauert  correction  (+0.00179)  (section  5.1)  has  been  applied  in  respect  of  tunnel 
constraint  effects,  giving  the  =  0.01069  total  entered  in  Table  3.  A  further 
correction,  for  strut  interference  on  the  Cp^  term,  is  also  appropriate  -  giving  an 
induced  drag  coefficient  Cp^  =  0.01079  .  (Table  3  also  gives  the  corresponding  data  for 
the  configurations  investigated  subsequently  and  reported  in  Refs  3  and  4.) 

b. b  Evaluation  of  results 

The  induced  drag  calculations  discussed  above  have  been  repeated  for  the  three 

1  3 

wing-alone  cases  investigated  at  75  m/s  (Table  I)  ,  and  the  results  are  plotted  in 

Figs  32  and  33.  Since  these  data  were  insufficiently  detailed  in  the  strut/wing  junction 

regions  no  strut  interference  Cdjj.  term  could  be  estimated,  and  the  value  for  the  wing 

plus  body  case  has  been  used  (Cdtt  =  -0.0001  at  Ct  =  0.491)  in  these  three  cases. 

.  lnt  c  2 

The  Cp  estimates  are  approximately  proportional  to  Ct  ,  and  the  induced  drag 
IVc  c 

values  are  consistent  with 


where  the  minimum  {ie  zero)  induced  drag  occurs  at  CT  =  -0.010  and  the  vortex  drag 

Lcv 

factor  k  =1.115. 
v 

On  both  the  above  figures  the  wing  plus  body  datapoint  is  closely  consistent  with 

the  lines  defined  by  the  three  wing  alone  datapoints.  Assuming  the  induced  drag  to  have 

14..  v 

a  parabolic  lift  dependence  ,  it  is  not  possible  to  deduce  both  k^  and  ^Lcv  on 

basis  of  a  single  traverse-derived  Cp^  value  (Cp  being  defined  by  a  mechanical  balance 

measurement),  for  the  wing  plus  body  case.  It  would  however  be  reasonable  to  assume  that 

V 

the  ^Lcv  vaiue  is  not  materially  altered  by  the  introduction  of  a  body,  and  accordingly 
estimate  k^  =  1.128  for  the  wing  plus  body  case.  This  result  is  consistent  with  the 
expectation,  from  standard  inviscid  theory,  that  the  introduction  of  a  body  modifies  the 
spanwise  lift  distribution  such  that  k^  increases.  It  can  be  assumed  that  the  25% 


Reynolds  number  disparity  between  the  wing-alone  and  the  wing  plus  body  cases  has  no 
'•iii'it  on  the  vortex  drag  characteristics  of  the  wing. 

niven  the  assumption  that  the  ininimum-vort  ex-drag  lift  for  the  wing  plus  body  case 
is  C[  =  -0.010  ,  and  the  additional  assumption  that  the  strut-on-model  interference 
on  the  total  drag  (not  included  in  the  mechanical  balance  data  shown  in  Fig  4 b )  is 
insignificantly  lift-dependent  in  the  range  0  <  Cp  <  0.65  ,  then  one  can  estimate  the 
litl  dependence  of  the  or  profile  drag  term.  The  use  of  Fig  4b  and  the  relations 
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produces  the  profile  drag  variation  shown  in  Fig  34,  where  Cn.  is  the  fully  corrected 

■  •  v  t  r  e  f  „ 

vertex  drag  at  the  balance-derived  minimum  total  drag  point  Cp  =  0.06,  =  0.0203), 

and  Cp  the  corresponding  profile  drag.  It  is  unlikely  that  the  assumptions  with 

1  re  f 

regard  to  the  negligible  lift  dependence  of  the  interference  effects  are  sufficiently 
invalid  for  the  basic  form  of  Fig  34  to  be  affected.  The  basic  conclusion  from  the 
Figure  is  that  the  minimum  profile  drag  for  the  wing-body  case  occurs  at  significantly 

.  V  ✓ 

higher  Iitt  (Cj^  =  0.35)  than  either  the  minimum  total  drag  condition  (C^c  =  0.06)  or 
the  minimum  vortex  drag  condition  (Cl  =  -0.01).  Over  a  limited  lift  range  the  profile 
drag  variation  can  be  very  approximately  represented  by: 


34  ♦  \(c  -  C  V 

CP I 


where  k^  =  0.4  .  As  for  the  parabolic  approximation  to  total  drag  variation,  a  wide 
range  of  k  values  might  be  chosen  according  as  to  the  closeness  of  approximation 
required  within  any  given  restricted  Cp^  range. 

Fig  35  shows  a  comparison  of  the  spanwise  distribution  of  local  profile  drag  con¬ 
tribution  (C$j(y))  measured  for  the  present  wing-body  case  (Up  =  60  m/s,  Cp^  =  0.49)  with 

that  for  a  broadly  comparable'^  wing-alone  case  (Un  =  75  m/s,  C,  =  0.48),  and  with  that 

1 3  C 

for  a  near  zero-lift  case  (Up  =  75  m/s,  Cp^  =  0.08).  In  all  cases  the  data  have  been 

scaled  on  local  wing  chord  (C^),  neglecting  the  small  spanwise  shift  between  the  point 

of  Cpj  measurement  on  the  traverse  plane  and  the  real  corresponding  -  ie  neglecting 

sidewash  effects.  In  the  outboard  region  ' B '  the  wing-body  and  wing-alone  results  at 

Cp^  —  0.5  agree  well;  the  rapid  fall  in  Cp^  close  to  the  tip  (more  than  proportionate 

to  the  planform  taper  effect  on  )  is  common  to  both  sets  of  results  and  indicates 

substantial  entrainment  of  wing  boundary  layer  and  wake  shear  layer  fluid  into  the  tip 

vortices.  In  region  'A',  for  the  two  Cp^  =0.5  cases,  the  wing-body  profile  drag 

contributions  tend  to  be  lower  than  those  for  the  wing-alone  case  -  the  disparity 

increasing  with  approach  to  the  spanwise  location  of  the  wing-root.  This  effect  can  be 

associated  with  the  strong  inwash  flows  generated  on  the  wing-bodv  conf igurat ion ,  by 

wing-root  vortex  effects  and  flow  convergence  over  the  rear  fuselage.  The  disparities 
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in  C]}j  distribution  between  the  two  Clc  =0.5  cases  in  the  mid-semispan  region  are 
dominantly  attributable  to  entirely  arbitrary  differences  in  the  location  of  the  wake 
analysis  regions  relative  to  the  main  features  of  the  wake  in  the  two  cases. 

In  the  near  zero-lift  case  the  outboard  of  the  two  maxima,  discernible  in 

the  other  two  cases  is  absent.  This  indicates  that  the  secondary  flow  normally  generated 
at  the  strut/wing  junction  (Fig  18)  is  also  absent  -  ie  that  the  effect  is  lift  dependent. 
The  low  levels  of  shown  by  these  data  indicate  a  considerable  net  increase  in 

profile  drag  for  the  wing-alone  configuration  between  Clc  =  0.08  and  CLc  =  0.48  . 

The  following  two  figures  extend  the  comparison  between  the  Ref  13  wing-alone 
results  and  those  from  the  current  wing-body  tests.  Fig  36  shows  that  there  is  a  con¬ 
sistent  trend  of  peak  with  lift,  the  peak  value  being  consistently  smaller  in  the 

port  vortex  core.  The  inverse  effect  is  apparent  in  the  results  (Fig  12c)  for  the 

wing-body  case,  and  it  may  be  recalled  from  the  table  in  section  5.2  that  in  this  case 
the  port  vortex  circulation  was  of  marginally  greater  magnitude  than  the  starboard. 

Fig  37  shows  that  the  vortex  core  span  for  the  wing-body  case  is  somewhat  less  than  for 
the  wingralone  configuration  at  the  same  lift,  though  clearly  part  of  the  difference 
relates  to  the  traverse  plane  being  1.6  tip  chordlengths  further  downstream  in  the  former 
case . 


7  EXPERIMENTAL  RESULTS  -  WAKE-TRAVERSE  LIFT 
7. 1  General  comments 

The  basic  procedure  for  calculating  the  model  lift  has  been  summarised  in  section 
3.3,  section  3.4  describing  the  evaluation  of  the  adjustment  for  the  strut  circulation 
effects  not  considered  in  the  original  analysis'.  This  adjustment  contributes  <.6%  of 
the  final  traverse-derived  lift  (Cj^ t  =  0.485,  uncorrected  for  tunnel  constraint  effects), 
which  falls  1.2%  short  of  the  corresponding  balance-derived  figure  (Table  2).  It  should 
be  noted  that  this  adjustment  does  not  imply  an  attempted  correction  for  strut  inter¬ 
ference  effects,  either  on  the  total  lift  or  its  spanwise  distribution. 


Fig  7  shows  the  spanwise  distribution  of  local  vertical  component  momentum  defect 
measured  in  the  traverse  plane.  In  global  terms  this  is  equivalent  to  the  configuration 
lift  but,  because  of  crossflow  effects,  the  distribution  shown  does  not  necessarily 

reflect  the  actual  spanwise  lift  distribution  on  the  wing  itself  -  particularly  in  the 
centresection  region.  The  total  lift  contribution  from  the  port  wing  wake  analysis 
regions  is  1.7%  greater  than  the  equivalent  starboard  contribution,  and  the  tip  vortex 
circulation  data  tabulated  in  section  5.2  shows  a  2.6%  discrepancy  of  consistent  sense 
in  vortex  strength. 


Fig  8  shows  how  the  local  strut  circulation,  induced  by  the  model,  varies  as  a 
function  of  distance  from  the  model  and  of  model  lift.  It  is  clear  that  the  distribution 
of  circulation  on  the  strut  is  strongly  influenced  both  by  the  strut  guard  leakage  plume 
and  by  the  strut/wing  junction  vortex.  It  is  also  clear  that  the  presence  or  absence  of 
a  fuselage  has  a  strong  effect  on  strut  circulation  distribution. 


27 


7 . 2  Evaluation  of  results 

A  comparison  of  balance-derived  and  traverse-derived  lift  data  (Table  2)  gives 
agreement  to  1.2%  in  the  present  case,  after  extraction  of  the  buoyancy  induced  lift  tares 
effect.  (As  with  the  buoyancy  induced  drag  tares  effect,  it  is  implicit  that  the  effect 
is  not  modified  by  introduction  of  a  model.)  A  similar  order  of  agreement  has  been 
obtained  for  the  Refs  3  and  4  results.  Given  the  uncertainties  involved  in  the  strut 
circulation  adjustment  and  the  imprecision  inherent  in  graphical  integration  techniques 
the  order  of  agreement  achieved  is  Considered  extremely  satisfactory.  On  the  basis  of 
these  three  cases  the  wake- trave rse  technique  seems  entirely  valid  in  lift  derivation 
terms. 

8  CONCLUDING  REMARKS 
8. 1  Experimental  results 

Detailed  traverse  measurements  have  been  made  in  the  wake  of  a  swept-wing  model  of 
wide-bodied  transport  type,  at  a  corrected  C^  of  0.491,  and  the  results  analysed  in  the 
context  of  the  Maskell  analysis  -  giving  the  drag  data  summarised  below. 


CD  *  '°4 

Net 

measured 

CD 

Auxiliary 

terms 

Vortex 

core 

adjustment 

Total 

tunnel 

corrections 

Correction  for 
support  rig 
interference 

Fully 
co  rrected 

CD 

Cdi 

228.7 

-0.5 

-4.4 

-12.8 

-4.4 

206.6 

cDn 

93.7 

-4.8 

+4.4 

+  13.6 

+  1.0 

107.9 

The  wake-traverse  data  has  been  compared  with  tunnel  mechanical  balance  data,  and  found 
to  be  1.2%  lower  on  lift  and  4.8%  higher  on  drag.  The  comparison  was  made  on  the  basis 
of  data  applicable  to  the  model  and  its  support  rig  combined,  before  application  of  any 
of  the  conventional  tunnel  corrections.  To  achieve  a  proper  basis  of  comparison  it  is 
necessary  to  remove  strut  guard  leakage  plume  effects  from  the  traverse  data,  and  to 
deduct  the  leakage  flow  buoyancy  effects  on  lift  (measured  as  AC^  =  0.006)  and  drag 
(estimated  as  AC^  =  0.00167)  from  the  balance  data.  It  is  also  necessary  to  correct 
the  traverse  data  for  the  effects  of  support  strut  circulation  not  allowed  for  in  Ref  1 

It  is  by  no  means  certain  that  the  effects  of  the  strut  guard  leakage  flows  have 
been  adequately  quantified  in  respect  of  drag,  and  that  the  4.8%  discrepancy  mentioned 
above  is  a  true  representation  of  the  accuracy  achievable  with  the  wake-traverse 
technique.  The  significance  of  the  leakage  flows  was  not  appreciated  during  the  main 
tests,  the  relative  positioning  of  the  guards  relative  to  the  struts  may  not  have  been 
satisfactorily  reproduced  in  the  supplementary  tests,  and  there  is  some  evidence  that 
the  buoyancy  drag  increments  are  sensitive  to  this  positioning. 

Consideration  of  the  and  Cdji  results  yielded  by  the  traverse  gives  no 

indication  that  they  should  not  be  literally  interpreted  as  'profile  drag'  and  'vortex 
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drag'  respectively.  However,  a  definitive  judgment  can  only  be  made  on  the  basis  of  an 

extensive  series  of  traverses  at  given  ,  with  variation  of  Reynolds  number  and 

traverse  plane  location.  This  should  certainly  be  done  on  a  'useful'  configuration  -  a 

swept-wing  -  without  the  complexities  of  high  lift  devices  or  their  deployment  mechanism 

I  3 

fairings.  Comparison  of  the  present  data  with  the  three  wing-alone  cases  suggests  that 
a  wing-body  configuration  is  no  more  difficult  to  analyse  than  a  wing-alone  configuration 
on  the  same  support  rig  (involving,  in  the  latter  case,  complex  centre  section  flows 
generated  by  the  necessary  sting).  Clearly,  the  resolution  of  the  current  uncertainty 
with  regard  to  the  strut  guard  leakage  flows  will  be  a  pre-requisite  for  the  above  series 
of  traverses. 

The  wing  plus  body  Cp^  data  is  consistent  with  a  vortex  drag  factor  kv  =  1.128 

(section  6.6),  assuming  a  similar  form  of  vortex  drag  variation  to  that  defined  by  the 

1 4 

three  wing-alone  cases.  That  is  with  a  parabolic  lift  dependence  -  and  a  near  zero-lift 

#  t  V 

minimum  vortex  drag  point,  at  ^Lcv  =  “0-01  in  fact.  The  above  kv  value  is  in  the 
expected  qualitative  relationship  to  that  defined  by  the  wing-alone  cases  (kv  =  1.115) 
where  the  spanwise  loading  will  be  more  nearly  elliptic.  However,  both  ky  values  are 
well  above  the  range  of  kv  variation  predicted  by  inviscid  linear  theory  (with  allowance 
for  wing  thickness  effects),  which  gives  ky  =  1.014  for  the  wing  plus  body  case  -  and 
at  most  k^  =  1.030  for  a  high-lift  configuration  (full-span  slats  and  40°  flaps).  The 
disparity  is  clearly  partly  attributable  to  the  neglect  of  viscosity  in  the  predictions, 
but  it  must  also  be  partly  attributable  to  the  assumption  of  a  planar  wake  in  such  pre¬ 
dictions.  (Even  at  the  traverse  plane  the  vortex  cores  already  lie  some  0.25  c  above  the 
centre  of  the  viscous  wake  in  the  mid  semi-span  region.)  The  precise  value  of  k^ 
inferred  from  cDjj  data  is  heavily  dependent  on  the  constraint  correction  applied.  The 
simple  Glauert  expression  has  been  considered  adequate  in  the  present  case,  and  its  use 
contributes  some  17%  to  the  final  fully  corrected  vortex  drag  value. 

The  above  k^  values  are  derived  assuming  that  the  'profile  drag'  measured  in  the 

vortex  core  (ACp  )  is  dominantly  a  consequence  of  the  dissipation  of  streamwise  vortic- 
rvc 

ity,  and  that  it  can  accordingly  be  transferred  to  the  Cp^j.  term.  If  this  had  not  been 

done  the  k^  values  for  the  wing  plus  body  and  wing-alone  cases  would  respectively  have 

been  1.095  and  1.082,  rather  nearer  the  theoretical  levels.  Clearly  the  validity  of  the 

ACp  transfer,  together  with  the  precise  constraint  correction  applied,  has  an  impor- 
Avc 

tant  bearing  on  the  interpretation  of  Cp^  measurements  and  will  require  detailed  study 
before  any  large  scale  wake-traverse  studies  are  undertaken. 

An  estimate  has  been  made  of  the  Cp^.  or  profile  drag  variation  with  lift  in  the 
vicinity  of  the  minimum  total  drag  point  (CLc  =  0.06),  on  the  basis  of  k  1.1 28  and  a 
minimum  vortex  drag  point  at  £r  =  ~®.0I  .  This  shows  a  minimum  profile  drag  lift  of 

V 

CpCp  =0.35  ,  though  it  is  unclear  as  to  the  extent  that  this  is  fixed  by  section  camber 
effects  and  the  extent  that  it  is  fixed  by  wing-body  junction  shear  flow  interaction 
effects. 

The  results  for  the  wing  plus  body  case  show  that  significant  mutual  interference 
effects  exist  between  the  model  and  the  main  struts.  Although  the  separation  of  the 
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struts-on-model  and  modeL-on-struts  components  of  the  net  interference  is  a  fairly  sub¬ 
jective  process  with  wake-traverse  data,  the  mechanism  of  the  interference  on  the  profile 
drag  is  clear.  A  secondary  flow  is  generated  at  each  outboard  strut/wing  intersection, 
presumably  because  the  strut  sectional  chordline  is  inclined  to  the  local  surface  velocity 
vector.  It  is  only  possible  to  identify  and  to  correct  for  the  first-order  interference 
effects  on  the  vortex  drag  term  -  in  terms  of  the  modification  of  the  spanwise  distribu¬ 
tion  of  shed  vorticity  produced  by  the  above  secondary  flow. 

The  present  experiment  is  insufficiently  elaborate  to  determine  the  pitch  strut 
interference  on  the  model,  though  it  is  clear  that  the  model  has  a  strong  adverse  inter¬ 
ference  on  the  pitch  strut  drag.  The  design  of  the  pitch  strut  top  and  the  strut/ 
fuselage  pivot/cleat  has  in  fact  been  shown  to  be  rather  poor. 

The  wake-traverse  measurements  reveal  two  pairs  of  wing/ fuselage  junction  vortices, 
the  upper  surface  wing-root  vortices  being  the  stronger.  The  circulation  of  these 
vortices  has  been  calculated,  together  with  those  of  the  vortex  pair  apparently  springing 
from  the  forward  fuselage  roof,  the  vortex  pair  generated  at  the  pitch  strut  top,  and  the 
two  strut/wing  secondary  flow  vortices.  The  three  pairs  of  fuselage  vortices  are  indi¬ 
vidually  quite  strong  (up  to  23%  of  the  tip  vortex  circulation)  but  make  small  individual 
contributions  (and  a  negligibly  small  net  overall  contribution)  to  the  or  ’vortex 

drag'  term.  Some  doubt  exists  as  to  the  precise  path  of  the  two  pairs  of  wing-root 
vortices  along  the  fuselage  sides. 

8. 2  Experimental  technique 

It  can  be  concluded  that  the  wake-traverse  drag  analysis  technique  has  considerable 
promise,  despite  the  9.8%  disparity  between  the  traverse-derived  and  balance-derived  total 
drag  figures.  The  evidence  is  that  this  disparity  originates  in  the  CDj  or  'profile 
drag'  term,  and  may  in  large  measure  be  attributable  to  uncertainties  in  establishing  the 
strut  leakage  flow  contribution  to  the  measured  mechanical  balance  drag. 

Given  the  reasonable  expectation  that  the  origin  of  the  above  drag  disparity  can 
eventually  be  firmly  identified,  one  then  must  consider  the  practicality  of  the  technique 
in  terms  of  the  resources  required  for  a  typical  drag  analysis  traverse.  It  is  unlikely 
that  the  requirement  for  measurements  at  5000-6000  points  can  be  reduced  significantly 
(although  experience  indicates  that  they  can  be  distributed  on  the  traverse  plane  to 
better  effect)  for  realistic  aircraft  configurations.  It  is  however  almost  certain  that 
large  reductions  could  be  effected  in  the  100+  hours  tunnel  running  time  currently 
required  for  a  traverse  -  by  improvements  in  the  traverse  control  software,  and  small  but 
significant  relaxations  of  the  yawmeter  nulling  criteria  which  need  not  lead  to  appreci¬ 
able  loss  in  the  accuracy  of  pitch  and  yaw  measurement.  Very  large  reductions  in  the 
data  analysis  workload  -  currently  some  1000  man-hours  per  traverse  -  can  be  obtained  by 
the  development  of  quite  simple  computer  software  for  the  automatic  plotting  of  processed 
data,  particularly  in  contour  map  form. 

In  short  the  wake-traverse  technique  for  drag  analysis  purposes  is  barely  practical 
as  implemented.  However,  some  relatively  modest  further  system  development  -  principally 
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in  the  computer  software  and  electronic  control  areas  -  the  technique  can  be  made  vastly 
more  cost-effective.  It  is  however  difficult  to  envisage  the  technique  ever  becoming  a 
practical  proposition  for  routine,  project-related  drag  analysis. 

In  assessing  the  usefulness  of  the  system  in  its  current  state  of  development,  one 
must  not  ignore  its  enormous  potential  for  general  flowfield  investigations  and  configur¬ 
ation  diagnostic  purposes.  These  areas  are  certainly  the  most  cost-effective  application 
of  the  system,  and  this  has  been  demonstrated  by  recent  vorticity  and  total-head-defect 
mapping  exercises  in  the  wakes  of  combat  aircraft  configurations  fitted  with  leading-edge 
strakes. 
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See  Tables  2  and  3  for  sunmary  of  lift  and  drag  resul 
See  Tables  4  and  5  for  detailed  traverse  drag  results 
Transition  fixed  on  upper  and  lower  surfaces 


LIFT  AND  DRAG  OF  MODEL  AND  ITS  SUPPORT  RIG  (UNCORRECTED) 
SUMMARY  OF  CASES  ANALYSED  IN  DEPTH 


<r 

m 

CO 

Oh 

x> 

<r 

8eap  asaaAea3 
-a^Bn  in  aoaaa  2 


8eap 

paAtaap  asaaAeaa 
-a^Bft  i«aoi 


Aaet  txne  iO0V 
OA3  343  }0  33N 


UU33  XjBtJTXnB 
IqD?  30  33N 


vC 

Oh 

m 

Oh 

O 

00 

CM 

vO 

O 

O 

O 

<3 

O 

Oh 

O 

Oh 

CN 

vO 

O 

O 

O 

vO 

fv 

O 

O 

O' 

m 

vO 

<r 

O 

d 

d 

CN 

00 

0 

Oh 

— 

— 

iA 

m 

V 

>h 

>S 

Ss 

T3 

T3 

TJ 

O 

O 

0 

.O 

X> 

X) 

T3 

XJ 

XJ 

m  C 

<y  c 

c 

vO  nj 

\0  CTJ 

—  00 

—  W> 

OO  GO 

C 

c 

VO  c 

(U  »H 

OJ  »H 

<f  *H 

<  3 

<  3 

3 

CQ  ^ 

cd  — * 

sc  ~ 

ajix  asaaAFaa 
ux  aoxaa  % 


109JJ3 

uot^FxnD:lTD  ^nj^s 
jog  pa^oaajo^ 


sajpq  3JTI 

SuiooBj^qns 


panoajaooufl 


pa^oajjooufl 


Mean  of  values  recorded  on  mechanical  balance  during  traverse 

Excluding  strut  guard  leakage  flow  effect 

See  Table  4(c),  totals  of  Cn  and  CnTT  contributions 
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Table  5 


ANALYSIS  OF  DRAG  CONTRIBUTION  OF  SUPPORT  RIG 
(a)  Traverse  measurements  in  starboard  strut  wake  (model-off) 


Source  of  C  contribution 

I 

CDX  X  '°4 

Strut  wake  region  (8)* 

Wing  wake  region  (4) 

15.3 

( i)  normally  exposed  strut  length 

5.2 

(ii)  strut  top  effect 

1.2 

Inferred  leakage  flow  effect 

6 . 7  t 

Mechanical  balance  Cn 

u  ta  re  s 

28. 4+ 

(b) 


Traverse  measurements  behind  complete 


support  rig  (model-on 


) 


t 


CD 

io4 

Source  of  drag  contribution 

port 

pitch 

starboard 

total 

strut 

strut 

strut 

Strut  wake  regions  (6,7,8)*^ 

16.5 

18.7 

15.9 

51 . 1 

Wing  wake  regions  (2,4) 

2.8 

- 

2.  7 

5.5 

** 

Fuselage  wake  region  (3) 

- 

8.8® 

- 

13.8 

Total  Cj)T  included  in  Cn 

1  utrav 

19.  3 

32.5 

*** 

' ■ ^*** 

3.3 

18.6 

70.4 

...  (strut  top  effect 

ln®  (leakage  flow  effect 

1.2 

6.7 

1.2 

6.7 

3.6 

16.7 

Subtracting  Cq 

Hnt 

-(-1.2)1 

-9. 8 1 

-  (-1 .9)1 

-6.7 

(model-on- struts) 

Mechanical  balance  C™ 

utares 

28.4 

27.2 

28.4 

84.0 

Notes : 

*  assuming  leakage  flow  gives  increment  only  for  guarded  portion  of  strut 

t  port  and  starboard  struts  assumed  to  give  identical  Cn 

6  utares 

**  strut  drag  increments  derived  by  extrapolation  of  curves  to  the  strut 

extremities 

***  assumed  values 

j-  deduced  as  the  residual  of  the  other  terms 
©,(D  see  Fig  1 3b(i) 
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A 

b 

b 

C 

c 


list  of  symbols 

wing  aspect  ratio 
wing  span 

effective  span  of  horseshoe  vortex  representing  wing: 

bl  =  L/(oU  r’) 
e 

cross-sectional  area  of  tunnel  working  section 
wing  mean  chord 


c 


T 


wing  tip  chord 


c(z) 


local  chord  of  main  strut 


S'  CL’  °M 

CD  *  CL  ’  CM 

c  c  c 


balance-derived  force  coefficients,  uncorrected 

balance-derived  force  coefficients,  corrected  for  all 
but  strut  interference  effects 


*  ^L 

c  cp  cv 


corrected  lift  coefficients  corresponding  to  minimum 
total,  profile  and  vortex  drags 


corrected  drag  coefficients  corresponding  to  above 


total  drag  and  lift  coefficients  from  traverse  data 
(uncorrected) 

corrected  induced  drag  coefficient  derived  from  traverse 
data 

profile  drag  coefficient,  local  contribution  per  unit 
spanwise  run  (mm-'),  local  contribution/ spanwise  mean 
contribution,  local  contribution  per  unit  area  (nnn-2)  - 
wake-traverse  derived 

vortex  drag  coefficient,  local  contribution  per  unit 
spanwise  run  (mm-'),  local  contribution/spanwise  mean 
contribution,  local  contribution  per  unit  area  (rran-2)  - 
wake-traverse  derived 

profile  drag  correction  for  blockage  effects 

'profile  drag'  increment  measured  for  vortex  core, 
accounted  to  Cp^  term 


ACDlIa’  ACDnb 


induced  upwash  auxiliary  terms  (Fig  bb),  drag 
coefficient  form 


CDt- 


Lmt 


CD 


Hint 


profile  and  vortex  drag  con t r ibut ions  from  support  rig 
interference,  in  coefficient  form 


°Di 


ref 


profile  drag  and  induced  drag  at  minimun  ’  ■  t a  1  drag 
condition 


local  lift  contribution  coefficient,  measured  in  wake- 
traverse  plane 


k 


constant  in  assumed  parabolic  law  of  total  drag 
variation 
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k 

P 


H 

H  . 
min 


AH 

P 

PO 

Pj>  P_*  P3*  P4 


P 


m 


6R 


5s 


U 


u 

e 

U* 


u  . 

circ 


u0 


v 


1  • 


o 


V 


2 


o 


1.1  Sf  OK  SYMBOLS  (continued) 

constant  in  assumed  parabolic  law  of  profile  drag 
variation 

constant  in  assumed  parabolic  law  of  induced  drag 
variation 

total  head 

minimum  total  head  in  wake  at  a  given  spuiwise 
posit  ion 

reference  total  head,  undisturbed  stream 
increment  of  total  pressure  defect 
static  pressure 

reference  static  pressure,  uruisturbed  stream 

'static  pressures'  measured  by  the  four  side  tubes  of 
a  five-hole  yawmeter 

tunnel  plenum  wall-static  pressure 

reference  dynamic  pressure,  undisturbed  stream 

dynamic  pressure,  local  dynamic  pressure  measured  by 
five-hole  yawmeter 

incremental  distance  in  anti-clockwise  sense  along 
liftloop  contour  (Fig  6c) 

incremental  cross-sectional  area  normal  to  tunnel 
stream 

streamwise  velocity  (x-axis  direction) 

reference  streamwise  velocity,  undisturbed  stream 

streamwise  velocity,  effective  value  produced  by  solid 
blockage  effects 

streamwise  velocily,  Betz  perturbation  velocity 
(JpU*2  =  Ho  -  (H  -  joU2) ) 

local  anti-clockwise  velocity  component,  tangential  to 
liftloop  contour  in  traverse  plane 

streamwise  velocity,  perturbation  contribution  produced 
by  solid  blockage  effects 

sidewash  velocity  and  sidewash  flow  angle  measured  in 
wake- traverse  plane 

downwash  velocity  and  downwash  flow  angle  measured  in 
wake- traverse  plane 

streamwise,  spanwise  (positive  to  starboard)  and 
vertical  (positive  upwards)  cartesian  co-ordinates. 
Origin  given  by  intersection  of  line  joining  T/E  tips 
with  vertical  C/^  plane 


x,  y ,  z 
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Ax 


a,  a 

’  t 

Yz 

Y 

r 


r 

Ar 


LK 


V  5. 


p 

6,  0 


B 


JL 


c/4 


^min’  ^max 


kB 


1> 


ra 

Y(y,  z) 


LIST  OF  SYMBOLS  (concluded) 

increment  of  streamwise  distance 

spanwise  position  of  main  strut  chordplane 

spanwise  position  of  intersection  of  the  wing  trailing- 
edge  with  the  fuselage 

spanwise  limit  of  application  of  strut  circulation 

correction  to  lift  data  (y'  =  y  ) 

s 

model  incidence,  five-hole  probe  pitch  error 

circulation  round  main  strut,  at  given  z 

circulation  round  arbitrary  contour  in  traverse  plane 

circulation  on  a  contour  defined  by  a  cut  through  the 
viscous  wake  and  the  liftloop  segment  on  its  right-hand 
side  (Fig  6c) 

maximum  value  of  on  a  semispan 

correction  for  strut  circulation  effects  (Fig  6c, 
section  3.3) 

denoting  an  incremental  quantity 

.  1 2 

constants  from  standard  tunnel  constraint  theory 
roll  tube  angular  position  error  for  five-hole  yawmeter 
yawmeter  pitch  angle,  pitch  angle  at  nulled  condition 
wing  sweep  at  quarter  chord 

streamwise  vorticity  measured  at  a  point  cn  the  wake- 
traverse  plane 

minimum  and  maximum  £  values  at  given  y 
fluid  density 

yawmeter  roll  angle,  roll  angle  at  nulled  condition 
five-hole  probe  yaw  error 

cross-flow  stream  function  on  traverse  plane, 
hypothetical  value  if  S(y,  z)  is  constant  to  infinity 
downstream 

cross-flow  stream  function  on  traverse  plane,  value 
actually  measured 
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Fig  2  Schematic  general  arrangement  of  model  mpport  rig 


Fig  5  Effect  of  transition  fixing  —  pitching  moments  about  mean  quarter -chord  axis 


Fuselage  section  at 
TE  wingroot  station 


i  showing  the  analysis  regions  (bounding  the  viscous  wake)  on  the  traverse  plane, 
locations  of  principal  model  features 


Fig  6b 
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Fig  6b  Summary  of  key  drag  expressions 


Fig  6c 


Liftloop,  following  edges 
of  analysis  regions 


'cire 
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V* 


13) 
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ArtK  =  strut  circulation  correction, 0  for  j  y/y'j  ^1 


Fig  6c  Division  of  traverse  plane  (x  =  const.)  into  analysis  regions, 
definition  of  liftloop  and  calculation  of  overall  lift 


Wing  alone,  75 m/sec 


Figs  9&  10 


Fig  9  Five-hole  yawmeter  probe  calibration,  relationship  between 
'measured'  and  'true'  dynamic  pressure 


Fig  10  Five-hole  yawmeter  probe  calibration,  relationship  between 
'measured'  and  'true'  total  head 
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Fig  13a  Distribution  of  CD|  contribution  (per  unit  run)  along  main  struts 


Fig  14  Variation  of  free -stream  conditions  across  the  empty  tunnel  working  section 


section  —  contours 
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Contours 


Fig  1 
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Fig  19  Locus  of  maximum  (negative)  observed  vorticity  and  its  value  -  port  viscous  wake 


Fuselage  contour 


Fig  22a&b  Measurements  in  fuselage  side  plane 


Fig  24a&b  Spanwise  distribution  of  downwasb  angle,  in  plane  of  tip-to-tip  line 
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Fig  25  Spanwise  distribution  of  downwash  angle,  below  the  viscous  wake  of  the  wing 


Spanwise  distribution  of  sidewash  angle,  in  plane  of  tip-to-tip  line 


Fig  27 


Fig  27  Spanwise  distribution  of  sidewash  angle,  below  the  viscous  wake  of  the  wing 


Fig  28a  Variation  of  streamwise  velocity  on  a  vertical  cut  and 
a  horizontal  cut  through  the  port  tip  vortex  core 


Fig  28b  Variation  of  sidewath  velocity  on  a  vertical  and 
horizontal  cut  through  tha  port  tip  vortex  core 


Variation  of  straamwise  vorticity  {£)  on  vertical  cuts  through  the  vortex  cores 
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Fig  29  Variation  of  downwash  angle  on  vertical  cuts  through  the  tip  vortex  cores 
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Fig  31  Contours  of  constant  total  head  defect  in  the  port  strut  guard  leakage  plume 
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Fig  32  Variation  of  ‘profile  drag'  contribution  measured  in  the 
vortex  cores  with  lift 


0.02 
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Fig  33  Variation  of  induced  drag  (corrected  for  constraint  and 
strut  interference)  with  lift 


Fig  35  Spanwise  distribution  of  profile  drag,  scaled  on  local  chord  (Cv) 
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